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ABSTRACT 


The presence of non-iinear behavior in potentiodynamic 
polarization plots has resulted in difficulty in determining 
the Tafel constants from such plots. A FORTRAN based program 
involving numerical differentiation techniques using a 
graphical display was used to determine the existence of the 
Tafel regions. 

Various alloys polarized in synthetic seawater and a 
5.57 NaCl solution were analyzed. Although severe 
concentration polarization often dominated the cathodic 
branches the techniques employed did allow for the selection 
of regions which approached linear behavior. The effects of 
concentration polarization in hindering the determination of 
Tafel constants were exemplified by the uncovering of a 
cathodic branch containing a small region where only 
actvivation polarization dominated followed by the onset and 
total domination of concentration polarization. 

A method of determining where the anodic and cathodic 
currents begin to dominate the potentiodynamic polarization 


Curve is introduced. 
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I. INTRODUCTION 


Recent studies at the Naval Postgraduate School have 
centered on the vibrational characteristics of certain high- 
damping alloys, [Refs. 1 - 8]. Most of these studies have 
been concerned with either the measurement of the damping 
capacities of these alloys and/or with the heat treatments 
necessary to attain the microstructures associated with the 
observed high-damping capacities. In addition, two separate 
studies were conducted in an attempt to determine the 
corrosion characteristics of these high damping alloys in 
marine environments. 

In the studies involving corrosion characteristics, both 
laboratory and actual marine environment exposure 
experiments were conducted. The actual marine environment 
exposure tests calculated corrosion rates based on direct 
weight loss methods. Тһе laboratory tests conducted were 
Linear Polarization and Potentiodynamic Polarization in 
synthetic seawater [Ref. 9] and a 3.5% NaCl solution [Ref. 
10] environments. Although the correlation between the 
laboratory and direct exposure methods was good, the 
individuals involved in both of these laboratory studies 
expressed concern with the difficulty in determining the 


Tafel constants necessary to determine corrosion rates. 


The actual experiments were conducted on the EG&G 
Princeton Applied Research Model 351 corrosion measurement 
system. The Model 351 utilizes a PAR Model 1000 
microcomputer featuring two Motorola 68000 microprocessors 
and a touch-screen input to create an environment which is 
well suited for conducting various corrosion tests. The 
System is, however, limited by the lack of a key board and 
more importantly by an operating system which is proprietary 
in nature. These two combine to make any analysis of 
experimental data, other than by the means proviced with the 
operating software, impossible unless the data can be 


transfered to another computer. 


11. DATA TRANSFER 


ест терта анс a ER ти з а ee es m 


The PAR 351 and its associated software has the ability 
both to plot the results of an experiment and to send just 
the data points to a printer via an RS 252 serial port. 
Although the software supports only Houston Instruments 
plotters, the printer option is configurable for numerous 
printers as long as the printer supports DC1 (XON, ASCII 17) 
and DCS (XOFF, ASCII 19) protocol. Other RS 232 
requirements such as baud rate, parity, word length, and 
stop bits are configurable as the user sees fit. 

During the initial data transfers, in the present work, 
Hayes Microcomputer Inc. SMARTCOM [I modem software was 
used to capture the data. This method was successful most 
of the time but since the two systems were not using modems 
but instead were directly connected via their serial ports 
some problems were experienced and editing of the data file 
after transfer was always necessary. The most serious 
detraction was the need to wade through various menus in 
order to accomplish the transfer. 

In view of the both the problems experienced with using 
SMARTCOM II and the neccessity for the user to be familar 


with the way SMARTCOM II works, it was decided to write a 


data transfer/capture program which would both speed 
execution and eliminate the need for user knowledge of 
communications software. 

Using the sample TTY data transfer program included in 
the manual for Microsoft GW-BASIC version 2 [Ref. 11] as a 
core a BASIC program was written to accomplish the transfer. 
The core was modified to eliminate the transfer portion and 
most importantly line feeds and carriage returns sent by the 
551 меге automatically deleted and a carriage return was 
inserted only at the end of a line containing data. This 
eliminated the problem of blank lines appearing within the 
data file corresponding to page headers and footers. Other 
modifications which were made involved checking for and 
correcting two formatting variations which resulted in 
compatibility problems between the transfered file and its 
use in FORTRAN application programs. The first involved a 
change from a floating format to a exponential format when 
the potential was equal to O.0000 volts. The 351 normally 
sends potential with 4 significant digits. When the 
potential is 0.0000 volts it is displayed as -ХХЕ-12. The 
formatting shift is not compatible with any FORTRAN formats 
as the total field width changes in addition to the minimum 
number of available digits. This occurance is Now trapped 
and the file modified to reflect the actual potential of 
0.0000 volts. The second formatting variation always 


occured and did so in а region which was of extreme interest 
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as follows. When the current density changes from negative 
to positive there is a range where the exponent is E-12. 
Unlike FORTRAN write statements which always reserve a space 
for the + or —- sign, ie. the + sign is not printed unless 
specified, the 351 used the space for the - sign if present 
or for the first digit if not. Again this resulted in a 
change in total field width at the occurance of -a.bcdE-12. 
The FORTRAN read statment would read this to be -a.bcdE-1. 
This was handled Бу simply deleting the least significant 
digit such that the number is read as -a.bcE-12 when this 
occurs. 

In 1ts final form the data transfer program requires the 
user to reply only a to a (R)eceive or (E)xit prompt and to 
enter a file name if the (Receive option is selected. 
Because the 351 can not write a selected configuration to 
disk the user is still required to properly configure the 


351 to the proper settings when powering up. These are: 


PROTOCOL S PRINT 
PARITY — NONE 
BAUD RATE - 2400 
STOP BIT ROME 
WORD LENGTH - 8 


Once these options are properly set on the 551 the 
desired data file should be loaded into the 351 and 


formatted into the linear option. At this point the user 
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selects the (R)eceive option on the Z248, presses the "PRINT 
DATA" prompt on the screen of the 351 and enters a filename 
on the Z248. 

Тпе final version of the data transfer program was 
compiled using Microsoft's BASIC compiler Version 1.0. By 
eliminating the writting of the data to the screen in 
addition to the file and using the compiled program vice the 
BASIC interpreter a file which will use OK of disk space 
can now be transferred in under 5 minutes as opposed to the 
20 minutes it took with ОМАНТСОМ ІІ. А listing of the 


program is included in Appendix B. 


12 


111. CORROSION RATES 


Нэн ——————————— CÓ (—ÁÓ— 8 8. 


а. | THEORETICAL EQUATIONS 

In general the corrosion behavior of a metal can be 
described by the use of the Nernst equation, activation 
potential, concentration potential and mixed potential 
theory. 

The Nernst equation is a modification of the Gibb's 
Free Energy equation and expresses the tendency of a 
corrosion reaction to proceed in terms of EMF based on the 


activities of the products and reactants of the 1/2 cell. 


E - E? — RT1og, (Q) 
nF 
Е - Electrode oxidiation potential 
(ø for reduction) 
E^ - Equilibrium oxidiation potential 


(9 for reduction) 
= Gas constant 
Temperature ^K 
Number of electrons taking part 
The Farday 
Activities of products/activities of 
reactants 


|| 


OQ TI3 — > 
|| 


The potential of the cell is the algebraic sum of 


Oxidation and reduction 1/2 cells. 


15 


Activation polarization describes the current density 


which exists when the cell is polarized. 


lanedic = leexp(+ZFaE/RT) 


le xc ss: = loGexp(-ZF(1-a)E/RT) 


lo = equilibrium exchange current density 
Z = Valance of the active species 
E = Polariztion about the equilibrium potential 


= fraction of E which applies to lancaic 
;1,F as above 


Da 


By manipulating the above equation the activation 


pou шр е can be written as: 


B logíio (15/12) 


7 
B 


2.SRT/aZF 


Although activation polarization occurs for both the 
anodic and cathodic reactions the cathodic reaction can also 
be effected by concentration polarization. When this occurs 
the reduction reaction occurs at such a rate as to decrease 
the concentration of ions in the solution near the surface 
of the metal. This results in a change in the potential of 
the 1/2 cell reaction since the effective activity in the 


Nernst equation changes. 
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The concentration polariztion с is expressed as: 


2 = = Du RT logio (1 = laz lu) 
nF 


operatinq current density 
limiting current density 


la 


The limiting current density is inversely proportional 


to the thickness of the diffusion layer. This thickness 
dependent on various geometric and environmental factors 
can often only be determined experimentally. The effect 
concentration polarization can be minimized by agitating 


solution, thus decreasing the layer thickness, or, as in 


15 


апа 


of 


the 


the 


case of oxygen reduction, by heating the solution to reduce 


the concentration of the dissolved gas. Both of these means 


of dealing with concentration polarization may significantly 


alter the desired corrosion environment and distort the 


correlation between the experimental results and the actual 


in-service corrosion rates. 


In the absence of concentration polarization a dynamic 


polarization diagram of a cell reaction would appear as in 


Figure 1. The slope of the anodic and cathodic curves 
de/d(logio ia) would be given by B as described in the 


activation polarization expression. 


Most cells are affected by concentration polarization 


such that the cathodic polarization = 277+ = = +7<- 


effect of this is shown in Figure 2. It is obvious from 


15 


The 


this that logiío (i14) approaches logis (ic) the linear 
behavior of the polarization curve is distorted and in cases 
where 1, and i are not sufficiently separate, a linear 
region may not exist. 

Mixed potential theory concerns concerns the immersion 
Of a single piece of metal into an electrolyte and the 
subsequent reactions. The theory consists of two simple 


hypotheses. 


1. Any electrochemical reaction can be diviced into two or 


more partial oxidiation and reduction reactions. 


2. There can be no net accumulation of electrical charge 


during an electrochemical reaction. [Ref. 12 : pp 314] 


It is these hypotheses that allows for the 
determination of two important parameters necessary for 
estimating corrosion rates. When the cathodic branch of the 
1/2 cell with the highest equilibrium potential intersects 
the anodic branch of the other 1/2 cell the corrosion 
potential Øceorr-r and the corrosion current density ico-- can 
be used to estimate corrosion rates. 

When the corrosion system consists of only a single 
metal in contact with a solution which contains a single 
reduction/oxidation reaction the overall reaction is the 
linear sum of the reduction currents and the oxidation 


currents. 
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If the solution contains more than one 
reduction/oxidation reaction the overall reaction may 
proceed as the linear combination provided that the 
equilibrium exchange current densities of the secondary 
reduction reactions are within one order of magnitude of the 
main reduction reaction. The secondary oxidiation reactions 
proceed in the same manner. As a result the secondary redox 
reactions affect either, both, or none of the branches of 
the overall reduction/oxidiation reaction. (Кет. 135: рр. 
199] 

If concentration polarization is not present both the 
anodic and cathodic polarization curves would consist of a 
of horizontial lines, representing the current added by the 
reaction at that potential, and slanted lines of some siope 
B representing the activation polarization of all the 
Oxidation or reduction reactions which are influential at 
that potential. In the presence of concentration 
polarization the possibility of introducing numerous поп- 
linearities into the cathodic branch is distinct. As the 
Overall reduction reaction is the sum of ali the influential 
reduction currents any region could consist of zones where 
the activation polarization of the ith reaction dominates 
the concentration polarization of the (1-1)th reaction. 

This can make the potentiodynamic curve extremely difficult 


to analyze. 
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B. DETERMINATION OF CORROSION RATES 

As previously described the presence of numerous 
electrodes and/or metals in the corrosion system can 
significantly alter the shape of the polarization curve. 
Despite this if it was possible to design a current 
measurement system that was capable of measuring only the 
anodic and cathodic currents the determination of icece- 
could still be simply made since 1l0gGi0 leo-re would be at the 
intersection of the anodic and cathodic branches of the 
polarization diagram. 

In the real world polarization curves take on a 'Y' 
shape as a result measuring instrument limitations and the 
complexity of the corrosion environment.  $-S-- can be 
measured directly since it is the equilibrium potential cf 
the system. As the potential is scanned towards Øeceorr the 
net current (iaás5oo:ic - lcatnoaic) 15 measured. Only when 
the potential is displaced far enough from 9$-4-- does one of 
the two currents overwhelm the other and does typical Tafel 
behavior take place. 

The presence of numerous 1/2 cell reactions, 
concentration polarization and the need for one of the 
currents to overwhelm the other can all combine to make the 
determination of leorr difficult. 

If concentration polarization is not present at some 659 
from é6-5-- the polarization diagram should adopt true linear 


behavior. The amount of 59 at which the curve should adopt 
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linear behavior is in itself open to disscussion. Ailor 
States that 59 should be at least 50mV from Øcorr LRef. 14: 
pp. 199]. Опсе the linear regions are found two methods can 
be employed to determine ic-o--. The simplest is to draw two 
straight lines from the linear regions towards eorr.» The 
two lines should intersect at ico--. A second method is to 
measure the slopes B (Tafel constants) of the linear regions 
of the anodic and cathodic branches. These values can be 
used in conjunction with data determined from a polarization 
resistance test to determine the value of i-o-- by use of 
the equation: 


icocec = В-ВМ 00) 


2.5(В.4В-)РН 
” 

If concentration polarization is predominant over the 
cathodic region of the polarization curve, the so-called 
‘knee’ method is often employed. In this method a tangent is 
drawn on the knee of the cathodic branch such that the 
angles formed on both sides of the tangent with the cathodic 
branch are equal. When this has been accomplished either of 
the methods discussed in the previous paragraph may be used. 

Many other methods have been developed to determine 
lco-- Or to determine corrosion rates directly from 
polarization data. Princeton Applied Research uses a Chi- 
squared minimization technique. This technique is included 


in the software package which operates their Model 551 


corrosion measurement system as PARCalc. By by their own 
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admission, since the program is terminated based on a change 
in the average value of the square of the relative 
deviations between two successive iterations, it is 
necessary to determine the Tafel region before starting the 
process. [Ref. 15: pp IX-10] 

LeRoy [Ref. 16:pp 10046-1012] postulated the 
determination of Tafel Slopes based on polarization 
resistance techniques. Without repeating the entire 
discussion it is sufficient to say his method of determining 
Tafel slopes based on expected values from activation and 
concentration polarization expressions, based on the 
assumption that the parameters necessary to evaluate those 
expressions were available, did not accurately describe the 
observed values. [Ref. 17: pp 1988-1989] 

It would also seem possible to solve the overvoltage 
expression for the terms in question. Recalling the 
overvoltage expression: 

$9 — $co-- = BelOGirolia/ia) + 2.3RT1IOGsao(l — 1а/1. ) 
nF 

The problem is that Be, 14 and ic all represent unknowns 
and the 1. term is located in such a manner as to complicate 


the solution. 
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The logíio (1 - 15/1,) term can be replaced by the series 


expansion: 


+ m А 


logio(1*X) = Z SL X — x + 
2 n 


(|х 
Ax 


where X = —-is/1l. 
The final equation to be solved would be: 
@ — eore = Bclogio(ia)-fi-logaio(is)*2.35(2(—-1)"**X^] 
n 
The difficulty with this approach is that a 
sufficient number of terms must be carried in the series in 
order to accurately represent the logio term. This causes 
two severe problems. The first is n*2 data points must be 
used to solve the system of equations. When one considers 
that ig is normally on the order of tA/cm=, exponentiating 
it to the nth degree can easily result in a matrix to be 
solved involving coefficients of SO or 60 orders of 
magnitude. The end result is that while the system can be 
solved such that the coefficients determined will result in 
the accurate representation of the actual overvoltage, the 
coefficients themselves have been influenced by round-off 
error and catastropic cancellation, and represent merely 
numbers. 
The other problem exists with attempting a solution of 
this form. As discussed previously with the possibility of 
mixed electrodes, the associated redox processes may or may 


not effect the overall behavior of the system. Even if the 
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limiting current densities i, were all known, it would still 
be necessary to determine their effect on the overall 
process in order to include them in the system of equations 
to be solved. As a result, for m different electrodes, m 
experiments involving just the electrode and the metal would 
have to be carried out just to determine their influence. 
Consider the use of such a method with synthetic seawater. 
At least 10 different chemical substances could effect the 
Overall reaction in addition to the major metal oxidiation 
and hydrogen reduction reactions. It should be obvious that 
if the limiting current densities are not known the size of 
the system of equations to be solved could easily exceed the 
number of available data points. More importantly, one 
would have to assess the advantages of performing a 
polarization in a mixed medium if the same test has to be 


performed separately on the individual components first. 


C. CALCULATION OF TAFEL CONSTANTS 

In order to determine the location, if one existed, of 
the anodic and cathodic Tafel regions two numerical 
differentation techniques were employed. These two methods 
were the Four Point Central Difference method and the use of 
a Cubic Spline interpolating polynomial. à graphical 


display of the derivatives allowed the user to determine the 
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region cf linearity. Once this region was determined a 
first order linear regression is performed on the original 
data within the region to determine the Tafel slopes. 

1. Four Point Central Difference Method 

When the data points are equally spaced the Central 

difference method may be used to approximate the derivatives 
at a point. The use of this method can result in 
considerable time savings over the cubic spiine method as 
the derivative can be calculated without having to resort to 
the solving of a system of equations for each point in 
question. The method is not as accurate as the cubic spline 
method and graphical representations of the derivatives will 
often contain more "noise" than those derived from the cubic 
spline method. 


The four point Central Difference formula for 


computing derivatives 1s: 


„ө = — шоолон RE o m mu 23 
d(logio Та) Gon — у) + (7—- — f2) 
Where: 
h = да - Фі-а 
fa = logaolia) at the Ith location from 


logiolla): 
The Central Difference method has an error term 
associated with it. In this case the error term is on the 
order of 1/h*. Although this may appear to be a substantial 


term, it should be remembered it is a possible error, the 
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method is only being used to determine the Tafel region and 
when used the Tafel constant is inverted so the error term 
would actually be h*. 

2. Cubic Spline Interpolating Polynomial 

When trying to numerically determine derivatives of 
experimental data which may or may not be equally spaced the 
use of interpolating polynomials is often suggusted. One of 
the most accurate methods is to pass a cubic spline thru the 
data points and then to differentiate the resulting equation 
to arrive at the derivative at a desired pcint. (Ref. 18: 
pp. 242] 

The advantage of using a cubic spline is that in 
deriving the system of equations to arrive at the constants 
a,b,c, and d which form the equation: 

y = ar(x — xr)> + Dr(x — xr)“ + cr(x — xr) FT RG; 

the system is formed in such a manner that the value 
of the function, first derivative and second derivative are 
the same for the pair of cubics which join at each point. 

This requires that: 

y° = Saxr(x — xy) + 261(х - жъ) + сү 
and уг = Gar(x — xx) + 25; 

In order to simplify the the mathematics involved 
the equations for the system are written in terms of the 


second derivatives of the interpolating cubic Sr- 
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The resulting interpolating cubic will then take on 


the form of: 


На-151-з1 “251 (Нг-з +5тг) + Н195т--1-Ф (ута Уг)- (УтгУх-а1) 
hz Нї-а 


where hx = x — Xr 


Solving the system leads to the followinq operations 


to determine the values of ag,Dxy,cix and di. 


ar = (бі-і - 5:)/6һһі 

бг = 61/2 

Сі - ((Уі-аі - У:)/Һһ:) - ((21һ:5і % һі5:..)/6) 
dx = У: 

The resulting system of n - 2 equations ап 5: 


involves n pairs of data points in order to generate the 


required number of equations. To arrive at the two 


additional equations to solve for S, and Snan, constraints are 


specified which pertain to the conditions at the ends of the 


curves. The three choices for the end conditions are: 


Sı = S8, = O, which implies that the end cubics 


approach linearity at their extremities. 


бі = Se, б, + 65,.-.і. This assumes that the end cubics 


approach parabolas at their extremities. 


5, 15 а linear extrapolation of S, and S2. S4 is a 


linear extrapolation of S,-2 and S4-;. 
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In this case since the interest is in determining 
the derivatives d¢@/d(logiao ia) pointwise, and not in 
developing an interpolating polynomial, it was decided to 
use the end conditions of S, = Sa = O, and to use enough 
data points as to force the inconsistences caused by the 
assumed end conditions away from the point in question. 

Four data points on each side of the point in question were 
used in determining the the derivative at the point in 
question. This resulted in the solving of a 7 x 7 system of 
equations for each data point. By locating the point in 
question in the center of the data set, the accuracy of the 
estimate involved was improved, and the mathematics involved 
in determining the derivative were simply solving for c;. 

4. Linear Regression 

Once the anodic and cathodic Tafel regions have been 
determined, the data points which comprise these regions are 
used ina least squares curve fitting routine to determine 


an equation of the form: 


ø = В(1одло (1а)) + 2, 


where B represents the Tafel slope. The values of a 
and b are found by solving the two following equations 
simultaneously. 

ВХ (109:.о(і-):) 2 + BEClOGraol(ia)1) - 2(104.о(14):) (9:) 

BZ(logao(l(ia):) + EN = 2(281) 


where N іс the number of points in the Tafel region. 


26 


When the PAR 351 Corrosion Measurement system scans 
the potential during a Potentiodynamic test the current is 
measured at intervals of 2mV or O.S5mV of potential, 
depending on the relative change in the current. The 
sampling rate is normally O.5mV when the current shows a 
substantial change and 2mV when it does not. This usually 
results in a sampling rate of O.5mV in the Tafel regions, 
although severe concentration polarization in the vicinity 
of izor- may result in a Limited number of points having 
O.SmV intervals on the cathodic branch. Two other aspects 
may also effect the potential intervals at which the current 
appears to have been sampled. 

When data is sent out the serial port, the 551 
normally uses a format similar to the FORTRAN Ga.b type. 
Not counting leading zeros, four significant figures are 
used. This results in a problem in the recorded potential 
values. When the potential is greater than or equal to 1.0 
volts, or less than or equal to -1.0 volts, only three 
decimal places are carried. às a result, when the values of 
$-o-- fall in this range and the scan rate is O.5mV, the 
recorded data seems to reflect that the potential has not 
changed for two successive current measurements. 

Since the potential is being scanned it is rather 
obvious when this occurs that the recorded data does not 


accurately represent the experimental data. Since this 
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pattern is easily recognized, the array containing the data 
can be manipulated tc more accurately reflect the 
experimental data. 

Occasionally the potentials recorded display no 
regular pattern and the data, although questionable in 
nature, must be used as recorded. 

When calculating the derivatives, the anodic and 
cathodic branches are performed separately. The 
calculations begin at the third data point away from EcS-- 
when the Central Difference method is used, and at the fifth 
when using the cubic spline method. The calculations 
continue as long as the data points continue to suggest a 
slope of of the expected sign (*/-). This automatically 
excludes the passivating region of the anodic section, as it 
should since it is not part of the Tafel region, and 
excludes the possibility of the program terminating 
execution if it were to calculate an infinite slope. In 
general this strategy worked as expected and provided enough 
data points to properly examine the Tafel regions. The only 
experiment in which it failed was for the SONOSTON sample 
tested in 3.5% NaCl solution. As can be seen by examining 
the Potentiodynamic curve of this sample (Figure 3), the 
cathodic branch exhibited behavior throughout which caused 
either method to terminate execution without generating any 


derivatives. The reason for this behavior is unknown. 


28 


IV. GENERAL OVERVIEW OF THE PROGRAM 


The program was developed so as to allow for its use 
without the supposition that the user was familar with the 
FORTRAN langauge or with personal computers in general. 
Numerous efforts were made to minimze the possibilty of the 
program prematurely terminating execution in the event of 
improper input or through floating-point operations 
resulting in invalid expressions. By taking advantage of 
the color capability of the monitor screen, colors are 
changed with respect to the type of input expected, and to 
display messages resulting from either input errors, data 
file format errors or a failure to properly follow the steps 
necessary for proper program execution. 

The program was written using Microsoft FORTRAN 
Optimizing Compiler version 4.0. Graphics were achieved 
using Microcompatibles, Inc.'s GRAFMATIC library. Plotter 
support was achieved using Microcompatibles, Inc.'s 
PLOTMATIC library. 

A summary of the main program and its subroutines is 
given below. The main program and its called subroutines 
are listed in Appendix C. The subroutines are presented in 


the order in which they appear in the main program. 
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а. START.FOR 
This is the main program. The name of the first data 
file to be used is entered. Following calls to DATAIN, 
CHECK and DATAM a menu consisting of & options is presented: 
1. Plot the potentiodynamic curve 
2. Generate Tafel Slopes using cubic spline 
5. Generate Tafel slopes using central difference 
4. Use another data file 
2. Overlay two potiodynamic curves 
5. EXTE 
Option i generates calls to GRAPH1 and PLOT 
Option 2 calls SLOPE1. Following this a first order 
linear regression of the anodic and cathodic branches using 
locations within the original array determined during a call 
То ОКАРН is performed and the results along with original 
potentiodynamic curve are passed to GRAPHI. 
Option 3 calls SLOPE. Upon return from SLOPE this 
option proceeds as Option 2 
Option 4 returns the user to the 3rd executable line of 
START. 
Option 3 calls DATAIN, CHECK and DATAM while retaining 


the data from the first file. 


Option & obviously exits the program. 


B. DATAIN.FOR 

This opens the input file and, using a formatted read 
Statement, generates the array containing potential and 
current density. Additionally this subroutine counts the 


total number of lines in the file. This is done to ensure 
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that formatting changes within the data file which were not 
corrected during the data file transfer will not have not 
caused the array generation to prematurely terminate. If a 
formatting error has been detected, an error message is 
displayed in red and the line number within the file is 
identified. This allows the user to exit the program and 


correct the error before continuing. 


C. DATAM.FOR 
This converts the absolute value of the current density 
to its logiíio value and determines the position of Eco-- with 


the array. 


D. CHECK.FOR 

This subroutine determines the number of data points in 
the anodic and cathodic branches to be used in the 
determining the derivatives when using the central 
difference method. This is done by checking log current 
density for values which might result in the calculation of 
an infinite slope (normally the onset of pitting or 
passivation) and for equal increments in the recorded 
potentials. Occasionally the recorded potentials will not 
reflect the polarization which has taken place. This occurs 
when the absolute value of the recorded potential is equal 


to or greater than 1.0. In these ranges the O.5mM voltage 
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increment falls outside of the total width of the field and 
is not recorded. As this is a recognizable pattern, at the 


users option the problem can be corrected. 


Е. 5 ОРЕЗБОН 
This is the first subroutine called when using the 
Central Difference method. It assigns graph titles, axes 


labels and calls DATADEL.FOR. 


F. DATADEL.FOR 
Using the points determined in CHECK, this calculates 
the numerical derivatives using the Central Difference 


method. 


G.  SLOUPEI.FOR 

This assigns graph titles and labels. It also 
determines the number of points, MI, to be used on both 
branches in the cubic spline method by checking only for 
values of logio current density which could result in the 
calculation of slopes of infinite value. It calls CSPLIN. 
separately for the anodic and cathodic branches. MI anodic 


is not necessarily equal to MI cathodic. 


H. CSPLIN.FOR 
For MI number of points, this generates the 7x7 matrix 
and the 7x1 vector of values described in the cubic spline 


section MI-9 times. As each matrix is formed it 1s solved 
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by a call to LINSY1.FOR. Upon return from LINSY1 the 7x1 
vector contains the solution which is used to calculate the 


derivative. 


I. LINSY1.FOR 
Using partial pivoting this subroutine solves the 
augmented 7x8 array by Gaussian Elimination. The original 


7x1 vector contains the solution [Ref. 19]. 


J. GRAPH1.FOR 

This is the graphics output of the selected main menu 
option. Through the use of the numeric keypad as softkeys 
the user can scale the axes as desired. The option to 
return to the original graph or the previous display always 
exists. When plotting the results of either numeric 
differentation method, only those points which consecutively 
fall within the desired region are plotted. Once the axes 
have been scaled such that graph contains only the linear 
region of the anodic or cathodic branches the position of 
the start and stop points within the original array 
containing the derivatives are saved. Since, when 
calculating the derivatives, the data point containing Ecorr 
is always used as a reference, ie. (depending on the method) 
the first derivative for each initial branch is either 5 or 
9 points away гот Ес, the start and stop points also 
represent the position of the Tafel regions within the 


original potential and log1O current density array. When 
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the graph is not of the numerical derivatives, any data 
point which falls within the desired axes is plotted. If 
desired GRAPHI calls PLOT.FOR for a hard copy output of what 


is currently displayed on the screen. 


K. PLOT.FOR 

This plots the graph displayed on the screen. Although 
many of the arguments passed from GRAPH1 are used directly 
as passed, this is not a ’screen dump’. The resolution of 
the pict: ғ” is much higher than that of the screen. AS a 
result what appears to be a one pixel jump on the screen 
will not plot as such. Two sets of tabulated data may be 
Plotted on each set of axes. At the users discretion 
different color pens may be used to plot the axes, first 


output, second output, and labels. 


L. CONFIGURATION 

In its final form, the program has been compilied and 
linked into an executable file which runs ona Zenith 248 
computer which is equipped with a 80286 microprocessor, an 
80287 Numeric Processor Extension, a EGA graphics card with 
64K video RAM, and a color monitor which supports the EGA 
color mode. The plotter used is an Houston Instruments DMP- 
40 series digital plotter. 

The plotter library is plotter series specific. Asa 
result of this, the use of a plotter other than a Houston 


Instrument's DMP-29 or greater would require obtaining 
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another library from Microcompatibles Inc. AS many of the 
subroutines contained within PLOTMATICS are similar to those 
called within GRAFMATICS, 1t can probably be safely assumed 
that little if any editing of the subroutine PLOT.FOR would 
be required. 

To use the program with a video configuration other than 


that previously described would require editing of the 


subroutines PLOT.FOR, GRAPH1.FOR, DATAIN.FOR, START.FOR and 
CHECK.FOR. DATAIN, START and CHECK would require editing 
only if the video configuration was monochrome. PLOT and 
СКАРН1 would require editing for any graphics confiquration 


other than EGA. 
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V. 


RESULTS AND DISCUSSION 


А. SYNTHETIC SEAWATER 


Table 1 below compares the Tafel constants determined by 


using the new program with those obtained by Escue [Ref. 


TABLE 1 


TAFEL CONSTANTS OBTAINED FROM 
SYNTHETIC SEAWATER EXPOSURE 


MATERIAL 


Ті-5ОМі 


1020 Steel 


3O4 Steel 


7075 А1 


Ее-Сг-Мо 


Ғе-Сг-А1 


Сч-Мп-А1-Ее-Ма 


Са-Мп-А1 


Cu-Zn-A1 


650 Bronze 


Ва 

0.1400 
0.0758 
0.0241 
0.3790 
0.1775 
0.1840 
0.0399 
0.0633 
0.0804 


0.0403 


€ From Escue [Ref. 9] 


Table 2 below compares the current densities and 


Вл,” 
0.174 
0.068 
0.178 
0.087 
0.121 
0.180 
0.048 
0.035 
0.062 


0.031 


Ве 

0.1155 
0.1481 
0.0207 
0.1584 
0.1226 
0.1369 
0.1579 
0.0652 
0.0785 


0.0724 


В-“ 
0.152 


0.114 


0.062 
0.141 
0.151 
0.120 
0.045 
0.075 


0.041 


9]. 


related corrosion rates of the samples tested in synthetic 


Seawater. 


current density as determined by the intersection of the 


Tafel lines. 


Tafel lines did not intersect at a single point and are 
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The first line for each material represents the 


In many of the cases the anodic and cathodic 


joined by a horizontal line. In these cases ico-- was 
determined by taking the average of the anodic and cathodic 
leorr at the point where the lines intersected a horizontal 
line form by @eorr- 


The second line represents the results obtained by 


combining the Tafel slopes with the polarization resistance 


data. 
TABLE 2 
CURRENT DENSITY AND CORROSION RATES 

MATERIAL 1с-2Т- ен CR (i. 

un / cm un / ст= mpy mpy 
Ті - 5ОУМі 1.245 1.754 0.57.5 0.799 
4.182 4.960 1.822 2.160 
1020 Steel 4.212 5.837 1.981 1.805 
20.325 17.201 9.580 8.110 
7075 A1 1.413 0.421 0.591 0.176 
3.067 0.994 1.281 O.415 
304 Steel 0.106 0.728 0.048 0.329 
0.157 1.190 0.071 0.537 
Fe-Cr-Mo 0.750 0.804 0.354 0.380 
2.294 2.060 1.085 O.974 
Ее-Сг-Ө1 О.822 О.ВЗО 0.362 0.365 
1.183 1.240 0.520 0.545 
Cu-Mn-Al-Fe-Ni 8.241 8.854 4.165 4.070 
9.662 10.400 4.882 5.260 
Са-Мп-А1 2.089 1.119 12151 0.595 
1.459 O.894 0.784 О.481с 
Cu-Zn-Al 25259 1.897 1.156 0.963 
5.780 3.230 1.919 1.640 
630 Bronze 1.762 1.442 0.901 O.757 
6.878 4.690 3.550 2.420 


* From Escue (Ref. 9] 
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1. Ti-50XNi (TINILG) 

These results compared favorably with those obtained 
by Escue. Figure 4 shows the derivatives as calculated by 
the cubic spline method. às can be seen in this figure no 
true linear regions are distinguishable in either of the two 
branches. When this occurs the "knee" method can be 
employed by examining the branches for a range where the 
derivatives oscillate about a linear midpoint. In this case 
both the anodic and cathodic branches showed this 
oscillation in the log current density range of —5.28 to - 
9.91. Figure 5 shows the resultant Tafel lines plotted on 
the original potentiodynamic curve. Although the curves 
intersect Ф.-с а «һа could be considered to be a single 
point the izor- determined in this manner differs greatly 
from that calculated by using the Tafel constants and the 
polarization resistance data. 

2. 1020 Carbon Steel (51020LG) 

Figures & and 7 show the results as obtained for the 
1020 carbon steel sample. Of all of the samples tested this 
one exhibited a cathodic Branch with the least tendency to 
behave in a linear fashion. If there was any linear 
behavior in the cathodic branch it occured at a logiío i-o-- 
Of about -5.2. Fhe anodic branch displayed a much longer 
lasting region around logis ic-o-- = — 4.5. This sample also 
showed the greatest disparity in the value of isore 


depending on the manner in which it was calculated. 
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3. 7075 Aluminum (AL7073LG) 

In this case the central difference method showed a 
linear region in both the anodic and cathodic branches. AS 
can be seen in Figure 8 the linear regions exist at the 
pronounced peaks in the derivatives. In this case the Tatel 
lines did intersect at a single point as seen in Figure 9. 

4. 304 Stainless Steel (SSTSO04LG) 

The calculated corrosion rates and the respective 
Tafel constants differed greatly from those obtained by 
Escue. Although the calculated corrosion rates correlated 
better to the actual weight loss methods, the need to 
terminate the calculating of derivatives before the onset of 
Passivity may have caused a linear region in the anodic 
branch to be omitted. These points can be seen by the sharp 
peak on the anodic branch on Figure 10, and the 
corresponding point on the potentiodynamic curve Figure il. 
The Tafel lines intersected at a single value of 10Gio leoer 
but this time there was a good correlation between the 
calculated and graphical ics--. 

9. Ғе-Сғ-Мо (VACRUILG) 

Again the presence of any truly linear regions is 
questionable, although both branches, Figure 12, show a 
range where the deviations from linear behavior are minimal. 
Figure 13 once again shows a theoretical intersection of the 


Tafel lines. 
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ө. Fe-tr-Al (VACRO2LG) 

This alloy showed behavior similar to the Fe-Cr-Mo 
alloy described above. What may be considered a truly 
linear region in the cathodic branch is shown in Figure 14 
just to the left of the sharp peak. The overlay of the 
Tafel lines on the potentiodynamic curve is shown in Figure 
15. The Tafel constants derived show a closer relationship 
to those calculated by Escue, than did a comparision of the 
Fe-Cr-Mo sample. 

7. Ма-А1-Ғе-Мі (50МО5ТІ б). 

SGNOSTON, Figures 16 and 17, displayed a linear 
anodic branch and, as did most of the alloys, a cathodic 
branch effected by concentration polarization. 

8. Си-Ма-бЯ1 (INCRLG) 

The result here again showed a linear anodic branch 
but in this case the concentration polarization did not 
overwhelm the cathodic branch in such a manner as to distort 
the region where activation polarization is predominant. 
Figure 18 shows the graph of all the derivatives and Figures 
19 and 20 show the graphs of the derivatives in what was 
decided to be the Tafel region. Figure 21 shows the 
original potentiodynamic and the Tafel lines. 

9. Cu-Zn-Al_ (DLALCLG) 

Experinced similar behavior as the Cu-Mn-Al sample. 

Once again a linear region was very apparent in the anodic 


branch. Figures 22 and 253. 
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10. 6390 Series Bronze (BRNZLG) 

This was the only sample of those tested in the 
synthetic seawater where the Tafel regions for the anodic 
and cathodic branches were found in the same region of log 
leor-. AS seen in Figure 24 concentration polarization 
again quickly dominates the cathodic branch. Figure 25 
displays the overlay. 

B. 3.5% NaCl SOLUTION 

In a Similar fashion Table 3 lists the Tafel constants 
and Table 4 the resultant current densities obtained from 
the samples tested in the 3.9% NaCl solution. In general 
the current results show a stronger correlation with the 
corrosion rates obtained by AGkthar [Ref. 10] than did the 
comparision of the current results with those obtained Бу 
Escue [Ref. 9] for the synthetic seawater. 

TABLE 5 


TAFEL CONSTANTS OBTAINED FROM 
5.5 4 NaC1 SOLUTION EXPOSURE 


MATERIAL B. В 2 < B< B< * 
304 Steel 0.1259 0.4117 0.0951 0.117 
7075 Al 0.0116 0.010 0.0585 2.56 
Fe-Cr-Mo 0.2821 0.254 0.1053 0.110 
Ее-Сг-А1 0.1203 0.214 0.1292 0.1285 
Cu-Mn-Al 0.0082 0.0137 0.0337 1.5797 
630 Bronze 0.0524 0.570 0.1472 0.321 


* From Akthar [Ref. 10] 
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TABLE 4 


CURRENT DENSITY AND CORROSION RATES 


MATERIAL beore паса CR ска 
un/cm-^ un /cm^ mpy mpy 
7075 Al 6.960 7.406 0.048 2.955 
10.460 11.700 1.842 4.685 
304 Steel 0.052 0.100 0.026 0.050 
0.012 0.020 0.006 0.010 
Fe-Cr-Mo 0.043 0.040 0.019 0.018 
0.609 0.613 0.274 0.276 
Fe-Cr-^1 0.493 0.460 0.234 О. 543 
0.916 1.180 0.428 0.551 
Cu-Mn-Al 2.864 6.062 1.454 3.078 
De 2 oe 10.800 2.666 3.483 
630 Bronze 0.093 0.129 0.048 0.066 
3.970 4.480 1.842 2.902 


ЖҒғотп Akthar (Ref. 10] 


1. 70795 Aluminum (967075) 
In this case the effect of concentration 
polarization made any estimate of the cathodic Tafel 
constant extremely difficult. As seen in Figures 26 and 27 
concentration polarization completely overwhelmed the 
cathodic branch. As a result only the "knee" method cuold 
be used, but its accuarcy is questionable. The anodic 


branch displayed typical linear behavior. 


The cathodic branch of this sample did display a 
linear region located between the sharp peaks in Figure 28. 
The anodic branch displayed numerous attempts to passivify 


as seen in Figure 29. 
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5. Ее-Сғ-Мо (УСМОРО) 

Similar behavior to the 3O4 Stainless steel above. 
The bandwidth exhibited by the cathodic branch near logio 
leorr Of about —- 6.5 Figure 30, typifies the knee in the 
curve. Again numerous attempts to passivify limited the 
range of data which could be use in the anodic branch Figure 
51. 

4. Fe-Cr-Al (VACAL2) 

Both the anodic and cathodic regions, Figure 32, 
displayed a region where the behavior could be considered to 
be linear. Although the Fe-CR-Mo and Fe-Cr-Al samples 
tested in synthetic seawater displayed similar behavior, the 
Fe-Cr-Al sample tested in the 3.5% NaCl solution did not 
exhibit the tendency to passivify, Figure 33, as did the Fe- 
Cr-Mo sample. 

9. Cu-Mn-Al (INCRMTE) 

This sample best typifies the problem associated 
with concentration polarization. On Figure 34 a linear 
region сап be seen in the cathodic branch at a 10Gio leorr 
of about -5.3. This region is followed by a rapid decrease 
in the slope of the curve indicative of concentration 
polarization. When the original potentiodynamic curve is 
viewed thru out its entire polarization range, Figure 55, 
very little can be determined due to the great changes in 


current density on the anodic branch. When the range in the 
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cathodic branch displaying the linear bevaior is scaled 
properly a textbook example of concentration polarization is 
seen, Figure 56. The overlay is shown in Figure 57. 
6. 630 Series Bronze (BRZ630PD) 

As seen on Figures 38 and 39 this sample experienced 
Similar behavior to the one above. All though the linear 
region directly before concentration polarization is not as 
well defined the cathodic behavior is dominated first by 
activation polarization and then by concentration 


polarization. 


C. COMPARISION OF CURRENT DENSITIES 

In many of the samples tested the corrosion rates as 
determined by the intersection of the Tafel lines were 
Significantly lower than those obtained by the use of the 
Tafel Slopes and linear polarization data. AS previously 
stated the Tafel regions should begin to dominant the 
polarization curve at an polarization increment of about 
SOmV from Ecor-r-- This implies that a linear region should 
not be present at increments less than this. With these 
tests a linear region was often found within this zone. 
Because of the .,eneral shape of a potentiodynamic 
polarization curve the resultant lines drawn from this 
region will always have a lower slope and thus intersect at 
a lower value of i. The disparity underscores the need to 


use the average of leorere as determined by the intersection 
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of the Tafel lines and that calculated by the Tafel slopes 
and linear polarization data when estimating corrosion 
rates. 

Another interesting feature was noticed when conducting 
some data manipulation of the potentiodynamic polarizations 
analyzed. When the absolute value of (9 —- 9$c-o--) was 
plotted against logiíio із іп many of the cases the anodic and 
cathodic branches began to deviate significantly at a log ia 
corresponding to the intersection of the two Tafel lines. 
This may be coinicidental but also may warrant further 
discussion. When a linear polarization resistance test is 
conducted the potential is scanned with a range of gdeorr %/- 
25 mV. The polarization resistance is determined by 
calculating the slope of the linear region about 1 = O 
A/cm2. If a linear polarization curve is plotted in the 
Same manner as above, while neglecting the singularity at i 
= 0 A/cm=, the deviation discussed would determine the end 
of the linear region. Returning to the equation used to 
determine icorre using the Tafel constants and the 
polarization resistance it can be shown that in order to the 
polarization resistance to be a constant that BaBe = 
CX(B.*tB-), where C is a constant of proportionality. Тһе 
current being measured during a potentiodynamic polarization 
is equal to the difference between the anodic and cathodic 
currents. Only when one of the currents begins to overwhelm 


the other would there be a substantial change in the 
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measured current with respect to a given overvoltage. This 
would correlate to the deviation between the anodic and 
cathodic branches ona plot of the absolute value of (96 - 
Фое Ме log (12e It also should correspond to the end 
of the linear region on the linear polarization curve since 
at this point the constant of proportionality would no 
longer be a constant as one of the two branches is beginning 
to overwhelm the other. It is postulated that a technique 
of this method could be used to determine the minumium 
overvoltage at which the search for a Tafel constant could 
begin. It 15 recognized that this does not account for the 
possibility of completely symetric anodic and cathodic 
branches but if this was the case there would be no problem 


in resorting to the +/- SOmV rule proposed by nilor. 
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VI. CONCLUSION AND RECOMMENDATIONS 


The FORTRAN based graphics program proved to be 
extremely valuable in determining the existence of linear 
Tafel regions particularly in those cases where 
concentration polarization dominated the cathodic branch. 

With the ability to transfer the experimental data to a 
programable computer the opportunity exists to further 
analyze the data. This could allow for the determination of 
which electrodes are effecting the main redox reaction when 
the environment consists of a solution of mixed electrodes. 
In order to accomplish this it is recommended that a single 
high damping alloy be subject to individual potentiodynamic 
polarizations containing a single component, at its normal 
concentration, of synthetic seawater. Through the use of 
graphical techniques it should be possible to isolate the 
influential redox reactions. 

In order to accomplish this it will be necessary to 
conduct all the experiments in solutions which contain the 
same amount, or at least a monitored amount, of dissolved 
Oxygen. This is one factor which was not considered in the 
previous two experiments and may have made the analysis more 


difficult than necessary. 
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In lieu of this approach an effort should be made to 
determine the limitinq current densities causing the severe 


concentration polarization. 
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APPENDIX B 
PAR 351 TO 2-248 DATA TRANSFER PROGRAM 


2-0 

SCREEN 

COLOR 4,1,1 

KEY ОЕЕ: 11 5:21 ОБЕ 

ГЕЕТМТ 6-7 

LOCATE 25,1 

PRINT STRING$(60," ") 

FALSE = O:TRUE=NOT FALSE 

XOFF$=CHRS (19) = XONS=CHRS (17) 

LOO LOCATE 25,20:PRINT "PAR 351 TO 2-248 DATA TRANSFER PROGRAM" 
110 LOCATE 3,20:PRINT "YOU MUST PROPERLY CONFIGUR THE PAR 351." 
120 LOCATE 5,20:PRINT “AFTER SETTING THE TIME AND DATE ON THE 351" 


355558535 


150 L(CAIE 7,20:PRINT "PRESS 'SYSTEM MONPCGER' " 
140 LOCATE 9,20:PRINT “PRESS CONFIAR SYSTEM" 
150 LOCATE 9,20:PRINI “SET “PROTOCOL” TO PRINT" 


160 LOCATE 10,20:PRINT “SET PARITY’ TO NONE" 

170 LOCATE 11,20:PRINT "SEI 'BAD RAE TO 2400" 

180 LOCATE 12,20:PRINT "SEI 'SICP BIT' TUO 1" 

190 LOCATE 13,20:PRINT "SET WORD LENGTH’ T0 8" 

200 COLOR 14,1,1 

210 LOCATE 15,20:PRINT “PRESS ANY KEY TO CONTINE" 

220 GOS = INEYS:IF GOt= "" GOTO 220 

250 TOES 

240 COLOR 4,1,1 
52025 


3 


INT “PAR 351 TO 2-248 DATA TRANGFER PROGRAM 


260 LOCATE 3,20:PRINT “PRESS ‘MAIN MENU’ ON TE 351" 

270 LOCATE 4,20:PRINT "PRESS 'RECALL/DISPLAY EXPERIMENT: " 

280 LOCATE 5,20:PRINT"PRESS 'DISPLAY EXPERIMENT'ON THE ‘COPY FROM DISK’ 
SECTION" 

290 LOCATE 5,20:PRINT “SELECT THE EXPERIMENT TO BE TRANSFERED" 

ЗОО LOCATE 7,20:PRINT "PRESS 'PLOT/FURMAT ' " 

310 LOCATE 8,20:PRINT “ENSURE THE X AXIS IS IN THE LINEAR SCALE" 

320 COLOR 2,1,1 

330 LOCATE 9,20:PRINT "IF NOT IN A LINEAR DISPLAY THEN PRESS" 

340 LOCATE 10,20:PRINT "'FOHMAT DISPLAY, ‘LINEAR’ , REVIEW DISPLAY’ , ‘PLOT 
FORMAT ' " 


350 COLOR 14,1,1 

360 LOCATE 14,20:PRINT “PRESS ANY KEY ТО СОМТІМ Е" 

370 GÜUS-INKEYS:IF GO$— "" GOTO 370 

зво а5 

390 COLOR 4,1,1 

400 LOCATE 25,20:PRINT “PAR 351 TO 2-248 DATA TRANSFER PROGRAM" 
410 SPEEDS < "2400" 

420 COMFILS ="COM2:"+SPEED$+" ,N,8,1,CS, DS" 

430 OPEN COMFILS АБ #1 


450 


470 


470 


510 


PRINT #1, XOFFS; 

(EN “SCORN: “ FOR ANPUT AS #2 

LOCATE 10,25:PRINT "(R)eceive a file OR (E)xit " 

LOCATE 11,37: INPUT TXRXS 

IF (TXRX$<>"R") AND (TXRX$<>"E") THEN 460 

IF TXRX$-"E" THEN 1:550 

ES 

COLOR 4,1,1 

LOCATE 25,20:PRINT "PAR 351 ТО 7-248 DATA TRANSFER PROGRAM" 
COLOR 15,1,1 

LOCATE 13,25:PRINT “PRESS ’PRINT DATA’ ON THE 551 THEN” 
LOCATE 15,30:PRINT “PRESS ANY KEY TO CONTINE" 

GOS = INKEYS: IF GOt="" GOTO 560 

(LS 

COLOR 4,1,1 

LOCATE 25,20:PRINT “PAR 351 TO Z-248 DATA TRAAGFER PROGRAM 
LOCATE 10,20:PRINT "THIS PROGRAM USES A RAM DISK TO SPEED 


EXECUTION." 


610 
620 
550 
640 
650 
660 
670 
680 
690 
700 
710 
720 
750 
740 
750 
760 
770 


COLOR 2,1,1 
LOCATE 13,20:PRINT “DO NOT ENTER THE DISK DRIVE OR THE FILE TYPE" 
LOCATE 14,20:PRINT "BUT ONLY THE NOME OF THE OUTPUT FILE" 

COLOR 14,1,1 

LOCATE 16,20:PRINT "ENTER THE OUTPUT FILE ME" 

LOCATE 18,25: INPUT FILS 

DSK$="D: " :DSKF IL$=DSK$4+F ILS 

OPEN "D:TMP" FOR OUTPUT 65 #3 

N = O 

as 

PRINT #1, ХОМ; 

IF LOC(1)-0 THEN COSB 820 

IF LOC(1)>128 THEN PAUSE=TRUE:PRINT #1,XOFFS 

P INPUTS (1,11) 

IF (AS=CHRS(10)) OR (AS=CHRS(13)) THEN ВЬ-СННЫ 15) :41-141 

IF (ӨВОСНБЫ 15) ) AND (A%<>CHRS(10)) THEN PRINT #3,A$; 

IF (J > 3) AND (B$-CHR$(13)) &OND (/&-CHR$(13)) THEN PRINT 


#5 „СННЪ( 14) 


780 
790 
800 
810 
820 
830 
840 
850 


IF (J > 3) AND (B$-CHR$(15)) АМО (б5-04НЫ (15) ) ТНЕМ Ве" "Е 
IF LOC(1)>0 THEN 730 

TF PAUSE THEN PAUGE=FALSE:PRINT #1, XONS; 
GOTO 720 

FUR І-і TO 5000 

IF t0C(1)<>0 THEN I = 9999 

NEXT I 

IF 1»9999 THEN RETURN 

OLOSE 85:01 S 

ЕФ="р: "каф 

OPEN "D: IMP" FOR INPUT 66 #4 

OPEN DSKF ILS FOR OUTPUT AS #5 

FOR I =1 TO N 

X$ = INPUTS$(44,#4) 

Y$ = INPUTS(2,#4) 


91 


950 
940 
950 


1000 
1010 
1020 
1050 
1040 
1050 
1060 
1070 
1080 
1090 
1100 
1110 
1120 
1120 
1140 
1150 
1160 
1170 
1180 
1190 
1200 
1210 
1270 
1250 
1240 
1250 
1260 
1270 
1280 
1270 
1500 
1510 
1529 
1540 
1540 
1550 


В% 
С% 
Dt 
ЕФ 
ЕФ 
G$ 
нъ 


MID5(X$,18,9) 

MIDS(B$,1,1) 

МІП%(В%,2,1) 

МІОФ(ВФ, 5,1) 

МІОФ(ВФ, 4,1) 

МІОФ(ВФ,5,1) 

MID$(3S,6,1) 

MID$(B$,7,1) 

J$ = MID&(BS,8, 1) 

K$ = MID$(B$,9,1) 

FX$ = "—0.000000" 

ТЕ (ЕФ = "E") THEN MID$(X$,18,9) 
ЈЕ (ЕФ = "E") THEN MID$(X$,18,9) 
IF (G&-"E") THEN MID$(X$, 18,9) 
IF (H$-"E") THEN MID$(X$, 18,9) 
IF (HS="E") THEN MID$(X5,18,9) — ЕХФ 

IF (I$-"E") THEN MID$(X$,18,9)-FX$ 

t^ (J&-"E") THEN MID5(X$,18,9)-FX$ 

I- (K$-"E") THEN MID$(X*,18,9)-FX$ 

F2$ = "E-12 " 

дъ - МІП%(Х%,42,1) 

21%-М106(Х%,43,1) 

82%-М10%(Х%,54,1) 

ASS = ASAIS 

IF (A3$ = "12") AND (42% - "-") TEN MID$(X$,39,5) = F2$ 
P$ = MIDS(X$,1,1) 

IF (P$ = "P") THEN FRINTSES, X 

NEXT I 

асеЕ #5:а5 

1-0 

COLOR 1,5,4 

BEEP 

LOCATE 12,30:PRINT "%82* TRANSFER COMPLETE e842" 

LOCATE 13,30:PRINI "PRESS ANY KEY TO CONTINE" 
GOS=INKEY$:1F БОФ="" БОТО 1270 

PRINT 81, ХОРЕФ 

CLOSE #1 

CLOSE #2 

CLOSE #3 

CLOSE #4 

COLOR 1,1,1 

GOTO 580 

END 


It 


H 
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APPENDIX C 
CORROSION PROGRAM 


THIS IS THE MAIN PROGRAM. WITHIN IT THE USER ENTERS THE FILE 
NAME WHICH HE/GHE WISHES TO ANALYZE. FOLLOWING THIS А МЕМ) 
APPEARS ALLOWING THE USER TO SELECT THE DESIRED OPTION. 


THE FOLLOWING VARIBLES ARE DEFINED 


СОАТА : THE 2 DIMENSIGNAL ARRAY FOR THE DATA FILE CONTAINING 
POTENTIAL AND CURRENT DENSITY 
AS ABOVE BUT ONLY USED WEN OVERLAYING TWO CURVES 
A TWO DIMENSIONAL ARRAY OF CDATA BUT NOW IN TERMS ОР 
POTENTIAL AND .0G10 QURRENT DENSITY 
XDAT =: A ONE DIMENSIGYAL ARRAY CONTAINING THE X VALLES 
(F THE FIRST CURVE TO HE GRAPHED OR PLOTTED 
YDAT =: A ONE DIMENSIONAL ARRAY CONTAINING THE Y VALLES 
CF THE FIRST CURVE TO BE GRAPHED OR PLOTTED 
XDAT1 : A ONE DIMENSIQWAL ARRAY CONTAINING THE X VALLES 
ОҒ THE SECOND CURVE TO BE GRAPHED OR PLOTTED 
YDATI : A СМЕ DIMENSIONAL ARRAY CONTAINING THE У VALLES 
OF THE SECOND CURVE TO BE GRAPHED OR PLOTTED 


CDATAL 
MDATA 


XI > THE SUM CF THE LOG OQURRENT DENSITIES. USED IN THE 
LINEAR REGRESSION FOR DETERMINING THE TAFEL CONSTANTS 
X2 > ПЕ ЯМ Ч THE SQUARES (F THE LOG CURRENT DENSITY 


INITIALLY USED IN THE SAME MANNER A XI BUT ONCE THE 
TAFEL. CONSTANT HAS HEEN DETERMINED USED AS THE LOG 
CURRENT DENSITY VALUE IN CREATING THE TAFEL LINES 


XY : ПЕ ЯМ G@ THE PRODUCTS F LOG QURRENT DENSITY AND 
POTENTIAL. USED IN THE SAME MANNER Аб XI AND X2 

Yi : THE SUM (F THE POTENTIAL. USED AS ХІ,Х2,ХҮ 

NS : THE NUMBER ОҒ POINTS USED TO GENERATE THE TAFEL SLOPE 

H2 > THE TAFEL SLOPE AG CA OE ATED IN THE REGRESSION 

B1 > THE INTERCEPT OF NE TAFEL LINE 

M > MENU RESPONSE 

I : DO LOOP COUNTER 

N > THE NUMBER (F POINTS IN THE DATA FILE BEING READ 

O1 > THE NUMBER OF POINTS IN THE FIRST DATA FILE, UNCHANGING 

м > THE NUMBER GF POINTS IN THE FIRST DATA FILE. USED IN 
THE SUBROUTINE GLOPEL AND CAN BE RETURNED WITH A 
DIFFERENT VALLE 

J = NÆ LOCATION ОҒ ECORR WITHIN MDATA 

41 > SAVE РБ J БІЛ РАБЗЕЮ АБ АЧ АНЗЕГЕМТ ТО Фа ОРЕ1 

02 > SAME AS 01 BUT USED UN Y WHEN UOVERLAYING DATA FILES 

J2 >: SAME AS J BUT USED ONLY WHEN OVERLAYING DATA FILES 

51 > THE POSITION WITHIN MDATA WHERE THE REQUIREMENTS CF 


THE CENTRAL DIFFERNCE METHOD FAIL. CATHODIC BRANCH 


< хх моми мии 345 я жж жж 6 жк »* € dX з » м 


Fer) г) СОЕ СЕО 


: СТЕ АБ S51 BUT FOR THE ANODIC BRANCH 

: ПЕ TOTAL NUMBER OF POINTS IN THE САТНООТС BRANCH 
WHIOH MEET THe REQUIREMENTS FOR USING THE CENTRAL 
DIFFERENCE METHOD 


54 >: SANE AS SS BUT FOR THE ANODIC BRANCH 

AL : THE POSITION WITHIN XDAT1 AND YDATI WHICH MARKS ІНЕ 
> BEGINNING OF THE TAFEL REGION 

a2 > SANE AS AL BUT THE END OF THE REGION 

AS > SERVES THE SANE PURPOSE 45 Al BUT WITH RESPECT ГО 
MDATA 

en > SANE (6 «s BUT AS А2 

CT > SAME СЄ А1 BUT WITHIN XDAT AND YCAT (CATHODIC) 

eZ > SAME 66 СЪ ВЛ THE END O THE REGION 

CS > SERVES THE SANE PURPOSE AS Cl BUT WITH RESPECT ТО 
MDATA 

C4 > SANE AS C3 АЛ AS C2 

NB : ЗЕТ ТО O. PASSED TO GRAPH AND PLOT ТО AVOID THE USE 
OF THE EMPTY ХПОТ1 A) YDATL ARRAYS WHEN THE ARE NOT 
USED 

NAME : THE NNE OG THE MAIN DATA FILE IN ШСЕ 

NAME? : THE NE OF THE SECOND FILE WHEN OVERLAYING 

TITLE2: THE SECOND HEADING FOR THE GRAPH AND PLOT. CHANGES 
DEPENDENT ON THE OPTION SELECTED 

TITLES: RESERVED FOR FUTURE USE IF DESIRED 

D1 : CHARACTER STRING USED ТО FORM TIILE2 

D2 : ІНЕ CHARCACTER REPRESENTATION OF B2 WHEN B2 REPRESENT 
THE CATHODIC TAFEL SLOPE 

DS > CHARACTER STRING USED TO FORM TITLE2 

D4 > SAME €% D2 BUT FOR THE ANODIC TAFEL SLOPE 

LIST = CHARACTER STRING USED FOR MESSAGES AND PROP IS 


REAL %4 XDAT (2000) , YDAT (2000) , XDAT1 (2000), 
+YDAT1 (2000) 

REALX8 MDATA(2000,2) ,MDATA1(2000,2) ,X1,X2,N5, XY, XI , Y1 
* ,CDATA(2000,2) ,CDATA1 (2000,2) 

ІМТЕСЕКЖ2 1,J,N,A1,A2,C1,C2,01,NI ,J1,S1 ,S2,S3,S4,02,J32,N8 
INTEGER#2 A3,A4,C3,C4 

CHARACTERSS МАМЕ NAVE, TITLE2Z&32, YNAYES5S,LIST 450, PTTLES*®50 
CHARACTERSB D2,D4,D1%5,D35%6 

ТЕРЕЛ =” 

УМЕМЕ = ' · 


THE SCREEN 1S QLEARED AND THE FILE MTE IS ENTERED. ПЕ USER 
CAN A S0 EXIT THE PROGRAM (R GET A LIST OF EXISTING DATA FILES. 
THIS PROMPT APPEARS IN PURPLE WHICH IS THE COLOR USED WHENEVER 
A CHARACTER RESPONSE IS EXPECTED TO A РНОМРТ. IF A FILE NAE IS 
ENTERED WHICH DOES NOT EXIST AN ERROR MESSAGE WILL APPEAR IN RED. 


CAL OSMODE(3) 
CALL GOLEAR(O,S) 
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L3 


МЯГИЙ МЭ ee 


DIS = о 
CALL. GOV ( 15, 12) 
WRITE(6,805) LIST 
САН ойм0м(25,12) 
LIST = ‘ENTER ТНЕ ЕНЕ МАМЕ TO BE LEED’ 
WRITE (6,805) LIST 
CALL GOMDV(20,11) 
LIST = “А BLANK TO EXIT, OR "DIR" FOR A DIRECTORY’ 
WRITE (6,805) LIST 
CAL GQCMN(CS6,10) 
READ(6,800) NAME 
CALL ОСМСМ(40,9) 
ЈЕ (МОРЕ „ЕД. ' ') THEN 
СОТО 900 
ENDIF 
IF (NOME .EQ. 'DIR') THEN 
CALL GOLEAR(O,5) 
CAL СОСМОМ(1,25) 
PAUSE ‘ENTER "DIR *. /W" OR A BLANK TO CONTINE’ 
CAL GER (25,6,1,1,79,0,5) 
GOTO 15 
ENDIF 
CALL DATAIN(CDATA,N,NOIE , IER) 
IF (IER .EB. 1) THEN 
GOTO 10 
ENDIF 
CALL DATAM(CDATA,MDATA,N,J) 
CALL CHECK (MDATA,N,J,S1,S2,53,S4) 
CAL GSMODECS) 
CALL GCLEAR(O,2) 


Ј1 = Ј 
М = N-1 
01 = М 
Ni = N 


THE MAIN MENU IS DISPLAYED IN GREEN. THE COLOR GREEN IS USED 
WHENEVER AN INTEGER IS EXPECTED AS THE RESPONSE ТО 2 РКОМРТ. 

IF THE INPUT IS NOT AN INTEGER GR A CHARACTER AN ERROR MESSAGE 
IS DISPLAYED IN RED. 


CALL QCMN (15,22) 

LIST 2 ' 

WRITE(6,805) LIST 

CAL GOMOV (15,18) 

LIST = 71. PLOT ORIGINAL РОТЕМТТООУМАМТС ORE’ 
WRITE(6,805) LIST 

CALL GOMOV(15,17) 

LIST = '2. GENERATE TAFEL SLOPES USING CUBIC SPLINE METHOD’ 
WRITE(6,805) LIST 

CAL 9СМОМ(15,16) 

LIST = '3. GENERATE TAFEL SLOPES USING CENTRAL DIFFERENCE METHOD’ 
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C T ГУС) 


АНТТЕ (6,905) LIST 
CALL GOMWV(15,15) 


LIST = °4. USE ANOTHER DATA FILE’ 
WRITE(6,805) LIST 
LIST = '5. OVERLAY TWO POTENTIGDYNAMIC CURVES’ 


CALL GOMOV(15,14) 
WRITE(6,805) LIST 
CALL GOMOV(15,13) 
LIST = “6. EXIT’ 
WRITE(6,805) LIST 
CALL GCMOV(29,12) 
LIST = ‘ENTER 1,2,3,4,5 OR 6’ 
WRITE(6,805) LIST 
CALL ОСМОМ (39,10) 
READ(5,810,ERR = 22,IOSTAT = ЈА) M 
CAL QCMM (39,9) 
IF ((J4 .NE. O) „ОР. (М „СТ. 6)) THEN 
CALL ОБМОРЕ (5) 
CALL GOLEAR(O,4) 
LIST = · · 
CALL ОСМОМ (20, 17) 
WRITE(6,805) LIST 
LIST = 'УОЈ DID NOT ENTER AN INTEGER OF’ 
CAL GCMOV(20,16) 
WRITE(6,805) LIST 
LIST = ‘VALLE 1,2,3,4 ŒR 5.’ 
CAL ОСМОМ(20,15) 
WRITE(6,805) LIST 
LIST = ‘PRESS ANY KEY TO CONTINE’ 
CAL QCMOV(20,13) 
WRITE (6,805) LIST 
СА ОСМОМ (39, 12) 
CALI. GINKEY(IS,IA) 
GOTO 20 
ENDIF 


OPTION 1 HAS BEEN SELECTED. THIS GENERATES A GRAPH AND OPTIONAL 
PLOT (F THE ORIGINAL POTENTIODYNANIC CURVE. 


IF (M .EG. 1) THEN 

DO 25 I - 1,01 
= MDATA(CI,2) 
= MDATA(I,1) 


TITLES-' 
CALL. GRAPHL(XDAT , YDAT ,01,J1, NOE , XDAT1 , YDAT1,NB, YNGE , TITLE2 


+ ,A1,A2,C1,C2, TITLES) 
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КӨЛЕМ ент) 


POI 49708) 


ern 


БОТО 20 
ENDIF 


OPTION 2 HAS BEEN SELECTED. THIS GENERATES THE DERIVATIVES OF 
THE POTENTIODYNAMIC CURVE USING THE CUBIC SPLINE METHOD. 


IF (M .ЕП. 2) THEN 
А1 
a2 
C1 
C2 
31 
Ni 


счаоооо 


111111 


ë 


а ПРЕ! (МАТА, № „71 „МАМЕ „А „А2 „С1 „С2) 


STILL WITHIN OPTION 2. THE PROGROM CHECKS ТО ЕМ ҒЕ THAT THE 
TAFEL REGIONS FOR BOTH THE АМООТС 6ND CATFODIC REGIONS НАМЕ 
BEEN SELECTED. 


IF ((C1 .GE. C2) .OR. (A1 .GE. А2)) THEN 
CAL ОБМО0Е (5) 

CALL QCLEAR(O,4) 

LIST = ' ' 

CAL QOMOV(15,20) 

WRITE(6,805) LIST 

LIST = ‘YOU DID NOT SELECT BOTH AN ANDDIC’ 

CALL 9СМОМ(15,19) 

WRITE (6,905) LIST 

LIST = ‘AND CATHODIC BRANCH OR ELSE THE POINTS’ 
CAL GCMNV(15,18) 

WRITE(6,805) LIST 

LIST = "МЕРЕ NUT PROPERLY TRAPPED. IN ORDER TO’ 
CALL ОСМОМ (15, 17) 

WRITE(6,805) LIST 

LIST = ‘PREVENT THE PROGRAM FROM CRASHING YOU’ 
CAL GCMV(15,16) 

МЕТТЕ (6,805) LIST 

LIST = "МЕТ START AGAIN.’ 

CALL QOM0V(15,15) 

WRITE(6,805) LIST 

LIST = ‘PRESS ANY КЕҮ ТО CONTINE’ 

CAL ОСМОМ 15,13) 

WRITE(6,905) LIST 

CAL. GINKEY(IS,IA) 

GOTO 20 

ENDIF 


ALL REQGUIRMENTS HAVE BEEN SATISFIED. LINEAR REGRESSION IS 
PERFORMED ON THE TAFEL REGIONS AND THE REGILTS ARE GRAPHED 


55 


45 


УМАМЕ = ‘VOLTS’ 

ОС 55 I = 1,01 
храт(т) = мрата(1,2) 
YDAT(L) = MDATA(I,1) 


CONTINUE 

X2 = 

XI = O 

XY = O 

NI = eBS(C2-C1) +1 
YI = O 

C3 = J-(44C1) 

C4 = J-(4+C2) 


DO 40 I + CA,C3 


XI = XI + MDATA(I,2) 
X2 — X2 + (MDATA(1,2)**2) 
XY = XY + (MDATA(I,2) aMDATA(I,1)) 
YI = YI+ MDATA( 1,1) 
CONT INUE 
NS = FLOAT(NI) 
Bl — ((XY£XI)-(YI£X2) )/ ( (XIKXI) -(NS8X2) ) 


B2 = ((XY8NS5)—-(YIXXI) )/((X2ZXNS)—-CXI&XI) ) 
Y1 = MDATA(J,1) 

X1 — ((Y1-B1)/B2) 

X2 = MDATA(C4,2) 

DX — (X2-X1)/(520.) 


00 45 I = 1,50 
YDAT1(1) — (B2XX2)«B1 
XDAT1(I) + х2 


X2 = X2— DX 
CUNTINLE 

B2 = -B2 

CALL CONV(D2,B2) 
X220 

XI = O 

XY = O 

NS = O 

YI = O 

АЗ - J+4+ñ1 
A4 = 144432 
DO SO I = 85,64 


XI - ХІ + мргла(1,2) 
X2 + Х2 + (МПАТА(1,2) #2) 
ХҮ = XY + (MDATA(I,2)3MDATA(L,1)) 
YI = YI+ MDATA(I,1) 
CONTINE 
NI = ABS(A3-A4)+1 
NS = FLOAT(NI) 


B1 — ((XYX3XI)-(YIX*XX2) )/ CCXIXXI ) -(N53X2) ) 
B2 — ((XYXNSG)-(YIXXI))/((X213NSj —(XIXXI)) 
Ү1 = MDATA(J,1) 


X1 = ((Y1-B1)/H2) 

X2 = MDATA(A4,2) 

DX - (X2-X1)/(50.) 

DO 95 I - 51,100 
YDAT1(I) = (B2*X1)+B1 
XDAT1(I) - X1 


il 


X1 — X1+ DX 
25 CONTINE 
05 =  Ва- 


CALL CONV(D4,B2) 

X2 = (XDAT1(SOO)+XDAT1(SOL) )/2 

TITLE2-D1/ /D2/ /D3/ /DA 

МВ = 100 

CALL СРАРН1 (ХОАТ „УБАТ ‚01,31 МАМЕ , XDAT 1, YDAT1,NB, YNEME , TITLE2 
+ ,A1,02,C1,C2, TITLES) 


THIS PART ALLOWS THE USER TO SELECT EITHER ANOTHER ANODIC OR 
CATHODIC TAFEL REGION CR TO REVIEW THE GRAPH WITHOUT HAVING 
TO DETERMINE BOTH THE TAFEL REGIONS AGAIN. 


gy o0000 


CAL GCLEAR(O,2) 
LIST=° ° 
CAL QON (20,14) 
WRITE (6,805) LIST 
LIST = ‘IF THE GREPH DID NOT SHOW A GOOD CORRELATION’ 
CAL ОСМОМ (20, 15) 
WRITE(6,805) LIST 
LIST = ‘YOU CAN REPEAT TE CALCULATIONS FOR THE’ 
CALL ОСМОМ(20,14) 
WRITE(6,805) LIST 
LIST = ‘BRANCH IN GLESTION' 
CALL @СМОМ(20, 15) 
WRITE(6,805) LIST 
LIST = '1. REPEAT CALOLLATIONS' 
CAIL GCMN (20,11) 
WRITE (6,905) LIST 
LIST = '2. EXIT THIS PORTION’ 
CAL OMN (20,9) 
WRITE (6,905) LIST 
LIST = ‘ENTER A 1 OR A 2’ 
CALL GCMN(20,7) 
WRITE(6,805) LIST 
CAL GMN (25,6) 
READ(6,B10,ERR=56) M 
ТЕ (М „ЕО. 1) THEN 
вото 50 
ENDIF 
IF (M .EQ. 2) THEN 
вото 20 
ELSE 
СОТО 56 


ENDIF 
ENDIF 


THE CENTRAL DIFFERNCE METHOD IS TO BE USED. THE OPTION PROCEEDS 
IN THE SAME MANNER AS OPTION 2 


IF (M .EQ. 3) THEN 


А + 0 

62 = 0 

C1 = O 

С? = О 

CALL SLOPE (MDATA,01,J1,NAME,A1,A2,C1,C2,51,S2,S53,S4) 


IF ((C1 .GE. C2) .CR. (A1 .GE. A2)) THEN 

СА GSMODE CS) 

CALL ӨЛЕАВ(0,4) 

LIST = “ ' 

СА: ОСМОМ 15,20) 

WRITE(6,805) LIST 

LIST = ‘YOU DID NOT SOLECT BOTH AN £CNODIC' 
CALL ССМОУ( 15,19) 

ИНТТЕ (6,905) LIST 

LIST = ‘AND CATHODIC BRANCH OR ELSE THE POINTS’ 
CAL ОСМОМ ( 15,18) 

WRITE(6,805) LIST 

LIST = ‘WERE NUT PROPERLY TRAPPED. IN ORDER TO’ 
CALL ОСМОМ (15,17) 

WRITE(5,805) LIST 

LIST = ‘PREVENT THE PROGRAM FROM CRASHING YOU’ 
CAL ОСМОМ ( 15,16) 

WRITE(6,805) LIST 

LIST = ‘MUST START AGAIN.’ 

CAL ОСМОМ ( 15,15) 

МЕТЕ (6,805) LIST 

LIST = ‘PRESS ANY KEY TO СОМТІМ Е: 

CAL ОСМОМ (15, 15) 

WRITE (6,905) LIST 

CALL QINKEY( 13,14) 

GOTO 20 
ENDIF 

рі='Вс= ” 

УМАМЕ = ‘VOLTS’ 
DO 60 I = 1,01 


XDAT(1) = MDATA(I,2) 
YDAT(1) = MDATA(I,1) 

CONTINE 

X220 

XI = O 

XY = O 

NI = ABS(C2-C1) +1 

YI = O 

C3 = J-(2+C1) 


100 


70 


100 


Са = 4-(24С2) 
DO 7O I = С4,С-55 


XI = XI + MDATACI,2) 
X2 — X2 * (MDOGTA(CI,2) £2) 
XY = XY + (MDATA(CI,2) &"DATA (CIS  1)) 
YI = YI+ MDATACI,1) 
CONT INE 
NS = FLOAT(NI) 
В1 - ((ХҮЖХІ)-(УІЖХ2))/((ХІЖХІ)-(М5%Х2)) 
H2 — ((XY&NS)-(YI&XI) )/ C(X23NS)—-CXI XXI) ) 
\1 = МАТА, т) 
X1 — ((Y1-B1)/B2) 
Х2 = MDATA(C4,2) 
DX = (X2-X1)/(S0.) 


DO 80 Г = 1,50 
YDAT1(I) — (B2XX2)4B1 
XDAT1(I) = X2 


х2 = X2- DX 
CONTINE 
В? = -B2 
CALL. CONV(D2,B2) 
D3 = Вас“ 
X220 
XI = O 
XV = O 
NS = O 
YI = O 
AS = J+2+ñn1 
да = JAHA 
DO 90 I = 63,04 
XI 2 XI * MDATA(CI,2) 
X2 = Х2 + (MDATA(I,2)**2) 
XY = XY + (MDATA(I,2)#MDATA(I,1)) 
YI = YI+ MDATA(I,1) 
CONTINE 
NI = ABS(A2-A1)+1 
NS = FLOAT(NI) 
Bl — ((XYXXI)-(YI€X2) )/ ( (XI£XI )-(NS£X2) ) 
B2 — ((XYXNG)-(YI£ZXI)) / ((X21tNS)-(XI XXI ) ) 
YL = MDATA(J, 1) 
ХІ = ((Y1-B1)/B2) 
X2 = MDATA(A4,2) 
DX = (X2-X1)/(50.) 


DO 100 I = 51,100 
YDATL(L) = (B2*X1)+B1 
XDAT1L(I) = X1 
X1 = X1+ DX 

CONT INLE 

CALL CONV(D4,B2) 

TITLE2=D1//D2//D3//D4 

NB = 100 
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105 


CALL. GRAPHI (XDAT , YDAT ,01 ,J1, NOE , XDAT1 , YDGT1, NB, YNEYVE , TITLE2 
РОДЕ а Се ТЕПЕ) 
CALL GOLEAR(O,2) 
Я. 
CALL ОСМОУ(20,16) 
ИБТТЕ (6,905) LIST 
LIST + "ТЕ ТНЕ ВВАРН 010 NOT SHOW A GOOD CORRELATION’ 
CAL GCMOV (20,15) 
WRITE(6,805) LIST 
LIST = ‘YOU CAN REPEAT THE CALCLLATIONS FOR THE’ 
CALL ОСМОУ (20, 14) 
WRITE(6,805) LIST 
LIST = ‘BRANCH IN GLESTION’ 
CALL ОСМЕМ(20,15) 
WRITE (&,805) LIST 
LIST = '1. REPEAT CALCLLATIONS’ 
CALL QCMDV (20,11) 
WRITE(6,805) LIST 
LIST = '2. EXIT THIS PORTION' 
CAL GCMOV(20,9) 
WRITE(6,805) LIST 
LIST = 'ENTER A 1 RA 2’ 
CAL 9СМОМ(20,7) 
WRITE(6,805) LIST 
CALL GCMOV(25,6) 
READ(5,810,ERR=58) M 
IF (M „ЕО. 1) THEN 
СОТО 57 
ENDIF 
ТЕ (М „ЕП. 2) THEN 
GOTO 20 
ELSE 
GOTO 58 
ENDIF 
ENDIF 


OPTION 4 RETUHNG THE USER TO THE BEGINNING OF THE PROGRAM 


IF (M .EG. 4) THEN 
GOTO 10 
ENDIF 


OPTICN 5 IS THE OVERLAY. ONLY THE. POTENTIODYNCOMIC СПБМЕ5 CAN 
HE GRAPHED OR PLOTTED. THE ORIGINAL FILE IS STILL AVAILABLE 
FOR FURTHER MANIPULATION 


IF (М .ЕО. 5) THEN 
ТІП ЕЗ-: 
CALL GO E£RR(O,5) 
LIST = ' 
CALL. GOV (15,20) 
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ИНТТЕ (6,805) LIST 
CALL OOMOV(25,12) 
LIST = ‘ENTER THE FILE NOME TO БЕ (СЕП: 
WRITE(6,805) LIST 
CALL GOMOV(20,11) 
LIST = ‘A BLANK TO EXIT, OR "DIR" FUR à DIRECTURY' 
WRITE (6,905) LIST 
CALL GOMOV(364, 10) 
БЕЛО( 6,800) МАМЕЗ? 
CAL. GOM: N(4O,9) 
IF (NAME2 „ЕО. · ') THEN 
GOTO 21 
ENDIF 
IF (N FE2 .ED. 'DIR') THEN 
CALL GOLEAR(O,5) 
CAL GOKN(1,23) 
PAUSE "ENTER "DIR t. /W" св а BLANK TO CONTINE’ 
CAL GUSCRL(23,4,1,1,79,0,5) 
GOTO 105 
ENDIF 
DO 750 I = 1,01 
XDAT(I) = MDATA(I,2) 
YDAT(I) = MDATA(I,1) 
CONTINE 
CALL DATAIN(CDATAL,02,NAMEZ2, IER) 
IF (IER .EH. 1) THEN 
GOTO 105 
ENDIF 
CALL DATAM(CDATAL ,MDATAL ,O2,J2) 
CALL CHECK (MDATAL ,O2,32,51,52,53,S4) 
DO 760 I-1,02 
ХРАТ1 (1) = МОАТА1 (1,2) 
УБАТ1(1) = МОАТАТ (Е, 1) 
CUNTINLE 
TITLE? = MERAY '"//NGME//^ '//NGE2 
YNE = ‘VOLTS’ 
NAME2=NAME 
МАМЕ = С 
CALL. GRCPHI (XDAT , YDAT ,01 ,J1 , NOE , XDAT1, YDAT1,02, YNGPE , TITLE, 
А1,82,С1,С2,ТГПЕЗ) 
МАМЕ = МАМЕ? 
CALL GCLEAR(O,2) 
GOTO 21 
ENDIF 
IF (M .EB. 6) THEN 
GOTO 900 
ELSE 
БОТО 20 
ENDIF 
FORMAT (АВ) 
FORMAT (A50, V) 


10.5 





THIS SUBROUTINE READG THE DATA FILE GENERTATED BY THE 
MODEL 351 WHICH HAS BEEN TRANGFERFED ONTO A DISK 


~ + 


THE ARRAY WHICH CONTAINING POTENTIAL AND 
CURRENT DENSITY 

ІНЕ М?Ж (F THE DATA FILE 

CHARACTER STRING PRESENT IN JHE DATA FILE 
CHARACTER STRING PRESENT IN THE DATA FILE 
NUMBER OF DATA POINTS IN THE FILE, 
DETERMINED DURING THE PROGRAM 


CDATA 


NAME 
C1 
ша 
N 


КР е СЕ СЕ DJ 


SUBROUTINE DATAIN (СОАТА „М „МАМЕ , IER) 
REA 48 CDATA(2000,2) 
INIEGER$2, N,I 
CHARACTER C1%17,C2*6,NOME48,LIST#50,C3450 
OPEN (UNIT = 10, ERR = 200,FILE = NOWE, STATUS = ‘GLD’) 
DO SO N=1,2000 
READ (10,900,ERR=60) C1,CDATA(N,1) ,C2,CDATA(N,2) 
CONTINUE 
REWIND( 10) 
DO 70 I = 1,2000 
READ(10,%,ERR = 100) 
70 CONTINUE 
100 IF (I .NE. N) THEN 
CALL GSMODE(3) 
CALL GOLEAR(O,4) 
LIST = '* 
CALL ОСМОМ ( 15,20) 
МЕТЕ (6,905) LIST 
LIST = "ам ERROR HAS BEEN ENCOLNTERED IN’ 
CAL ОСМОМ ( 15,19) 
WRITE(6,905) LIST 
LIST = ‘THE FORMAT OF THE DATA FILE.’ 
CAL ОСМОМ ( 15,18) 
WRITE(6,905) LIST 
LIST = ‘THIS USUALLY INVOLVES A SWITCH BETWEEN’ 
CALL ОСМОМ ( 15, 17) 
МКІТЕ (6,905) LIST 
LIST = ‘FLOATING POINT AND EXPONENTIAL NOTATION’ 
CAL ОСМОМ ( 15,14) 
WRITE(6,905) LIST 
LIST = ‘AROUND ECDCR. ' 
CALL ОСМОМ ( 15,15) 
WRITE(6,905) LIST 
LIST = ‘TO CORRECT THIS PROBLEM EXIT THE PROGRAM’ 
CAL ОСМОМ (15,13) 
WRITE(6,905) LIST 
LIST = ‘AND USE THE DOS "EDLIN" FEATURE. TE’ 
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CALL GOMOV(15, 12) 
WRITE (6,905) LIST 
LIST = ‘PROBLEM GEGINS AT LINE #:' 
CAL OOMOV(15, 11) 
WRI TE(6,905) LIST 
CALL GOMOV(38, 11) 
WRITE(6,915) N 
LIST = ‘PRESS ANY KEY TO CONTINE’ 
CALL GOMOV(15,8) 
WRITE(S,905) LIST 
CALL. GINKEY ( I3, I4) 
GOTO 990 
ENDIF 
СОТО 990 
CAL QS"DDE() 
CAL GOLEAR(O,4) 
LIST = ° 
CAL. QCWN 1175,20) 
WRITE(6,905) LIST 
LIST = ‘THE FILE NAE ENTERED DXES NOT EXIST.' 
CALL ОСМОМ (15,19) 
WRITE(5,905) LIST 
LIST = "TO EXAMINE THE DIRECTORY OF POSSIBLE FILES’ 
CAL QOWN(15,18) 
МЕТЕ (6,905) LIST 
LIST = 'ENTER "DIR" WEN ASKED TO ENTER THE FILE’ 
CAL ОСМОМ ( 15, 17) 
WRITE(6,905) LIST 
LIST = 'NaME FUüLLOWED BY THE DOS COMAND.’ 
CAL ОСМОМ (6,905) 
WRITE(5,905) LIST 
LIST = '"DIR *. /W'' 
CAL GQOWM(CSS, 16) 
WRITE(6,905) LIST 
LIST = ‘PRESS ANY KEY TO CONTINE’ 
CAL GOMN(15,14) 
WRITE (6,905) LIST 
IER = 1 
CALL GINKEY(13,14) 
GOTO 990 
FORMAT (A17,G8.4,A7,G10.4) 
FORMAT (ASO, \) 
FORMAT ( I4, ) 
RETURN 
END 
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THIS SUBROUTINE MANIPLLATES THE DATA IN THE ARRAY CDATA 


THE PROPER FORM AND IDTENIFIES ECORR 


CDATA = THE ARRAY CONTAINING POTENTAIL AND CURRENT 
DENSITY 


MDATA - СТА CONVERTED TO POTENTIAL AND LOG 
CURRENT DENSITY 

N = NUMBER ОҒ DATA POINTS IN THE DATA FILE 

J = LOCATION WITHIN PDATA WHERE ECORR IS LOCATED 

I = DO LOOP COUNTER 

X1 = ABSOLUTE VALLE OF CURRENT DENSITY AT I 
USED FOR DETERMINING ECCRR 

X2 = ABSOLUTE VALUE OF CURRENT DENSITY AT I-1 
USED FOR DETERMINING ECORR 

X3 < 10610(Х1) 


SUBROUTINE DATAM (CDATA,MDATA,N,J) 
БЕЛ ЖӨ MDATA( 2000, 2) ,CDATA(Z000,2) 
INTEGER$2 N,J 

Х2 - 100000. 


CONVERT CURRENT DENSITY INTO LOG CURRENT DENSITY 
DO 100 I-1,N-2 


Xi = ABS(CDATA(I,2)) 
X3 = LOGIO(X1) 


IF/THEN TEST USED IN DETERMINING LOCATION OF ECOOR 


IF ( X1 .LT. X2) THEN 


1-1 
X2 = Xi 
ELSE 
X2 = X2 
ENDIF 
MDATA(I,1) = CDATA(I,1) 
MDATA(I,2) — X3 
CONTINE 
RETURN 


END 
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THIS SUBHÜUTINE PEHFOHMÉ TWO FUNCTIONS. ITS FIRST FUNCTION 
IS TO CHECK TO ENGLFE THAT THE POTENTIAL VALLES HFECDRDED 
АСО НАТЕ У КЕР ЕСТ THE ACTUAL POLARIZATION. THE SECOND IS 
TO DETERMINE THE NUMBER OF POINTS IN THE ANODIC AND CATODIC 
BRANCHES WHICH MAY BE USED IN THE CENTRAL DIFFERNCE METHOD. 


мрата : THAT ARRAY CONTAINING POTENTIAL AND LOG CURFENT 


DENSITY 

NS > THE NUMBER O POINTS IN THE CATHODIC ВЕРАСН WHICH 
CAN BE USED IN THE CENTRAL DIFFERENCE METHOD 

ма : THE NUMBER OF POINTS IN THE ANODIC BRANCH WHICH 


CAN ВЕ USED IN THE CENTRAL DIFFERNCE METHOD 
OTHER VARIBLES ARE USED MERELY IN CHECKING REQUIRED CONDITIONS 


SUBROUTINE CHECK (MDATA,N,J,N1,N2,N3,N4) 
REAL*8 MDATA(2000,2) 

INTEGER$2 1,J,N,N1,N2,N5,NA 
CHARACTERSSO LIST 


THE FIRST INCREMENT IN POTENTIAL FUR ME САТНООЈС ВЕРМСН 15 
DETERMINED 


Y = MDATA(J,1) 
Ү1 = MDATA(J-1,1) 

11 = ABS(ABS(Y1)—-ABS(Y)) 
71 = ANINT(Z1/0.0001) 


SUBSEULENT POTENTIAL INCREFMENTS, WHIUH M 5T EE EG XX. TO THE 
FIRST FFE CHECKED. LOG CURRENT DENSITY VALES АНЕ (4 50 СНЕСКЕО 
TO ENSURE THAT THEY WILL NOT RESULT IN THE CAL(ALATION (F ñ 
SLOPE OF INFINITE VALLE. 


DO 20 I = J,2,-1 
Y = MDATA(I,1) 
Yi = MDATA(I-1,1) 
X = MDATA(I,2) 
X1 = MDATA(I-1,2) 
Z = ABS(ABS(Y1L)-ABS(Y)) 
Z = ANINT(Z/0.0001) 


ONE OF THE TWO REGUIRED CONDITIONS HAS FAILED. 
IF ((Z .NE. Z1) .CR. ( X1 .LE. X )) GUTO 25 


CONTINE 
Ni = I+2 
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THE FIRST INCREMENT IN POTENTIAL FOR THE ANODIC BRANGH IS 
DETERMINED 


MDATA(J+1,1) 
MDATA(J+2,1) 

1 = ABS(X-Y) 

1 = ANINT(Z1/0.0001) 


X 
T 
Z 
7 


GUBSEQLENT POTENTIAL INCREMENTS, WHICH MUST BE EGUA ТО ПЕ 
FIRST ARE CHECKED. LOG CURRENT DENSITY VALUES ARE ALSO CHECKED 
TO ENSURE THAT THEY WILL NOT RESULT IN THE CALCULATION (F A 

Ч ОРЕ OF INFINITE VALLE. 


DO ZO I = J42,N-1 
Y = MDATA(I,1) 
Yl = MDATA(I+1,1) 
X = MDATA(I,2) 
Х1 = MDATA(IT1,2) 
Z = ABS(Y1-Y) 
Z = ANINT(Z/0.0001) 


ONE OF THE TWO REGUIPFED CUNDITIONGS HAS FAILED. 


IF ((2 „МЕ. Z1) .OR. ( X1 .LE. X )) СОТО 35 


CONT INGE 

N2 = I — 2 
N3 = J — N1 
м = ме - 9 


THE NUMBER OF POINTS AVAILABLE FOR EACH RANCH IS СОМРАНЕО 
WITH A MINUMILUM NUMBER OF 15. IF NETTHER HRANGH HAS ENOUGH 
AVAILABLE POINTS VARIOUS MESSAGES ARE DISPLAYED LISTING 
POSSTBLE CAUSES AND OPTIONS. 


IF ((NS АЕ. 15) „СВ. (NA .LE. 15)) THEN 
CALL GOLEAR(O,4) 
LIST = · 
CALL GOMOV(20,20) 
WRITE (6,805) LIST 
LIST = 'THE CENTRAL DIFFERNCE METHOD REQUIRES THAT’ 
CAL GOM0V(20,19) 
WRITE(6,805) LIST 
LIST = ‘THE CURRENT НЕ EARED АТ EQUAL DIFFERENTIALS’ 
CALL GOMOV(20,18) 
МЕТТЕ (6,805) LIST 
LIST = ‘OF THE APPLIED VOLTAGE. THE PAR MODEL 351’ 
CALL GOMOV(20,17) 
WRITE(6,805) LIST 
LIST = ‘RECORDS THE CURRENT AT INCREMENTS OF 2.0 mV’ 
CALL GOMOV(20, 16) 
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WRITE(&,805) LIST 
LIST = “UNTIL THE CHANGE IN CURRENT EXCEEDS A PRESCRIBED’ 
CALL ОСМОМ (20, 15) 

WRITE(6,805) LIST 

LIST = ‘LEVEL. AT THAT POINT THE 351 BEGINS RECORDING’ 
CAL GCMOV(20,14) 

WRITE (6,805) LIST 

LIST = ‘THE CURRENT AT INCREMENTS OF 0.5 ам. АТ ӨМ 
CAL GCMOV(20,13) 

WRITE (6,905) LIST 

LIST = ‘HIGHER APPLIED VOLTAGE THE CHENGE IN CURRENT’ 
CALL QCMOV(20,12) 

WRITE(6,805) LIST 

LIST = “DECREASES TO A POINT WHERE THE CURRENT IS CGAIN’ 
CALI. GCHOV(20,11) 

WRITE(5,905) LIST 

LIST = 'MEAARED АТ INCREMENTS OF THE £CPPLIED МО ТОБЕ 
CALL GCMOV(20,10) 

WRITE(6,805) LIST 

LIST = ‘EQUAL TO 2.0 aw’ 

CALL ОСМОМ(20,9) 

WRITE(6,805) LIST 

LIST = ‘PRESS ANY KEY TO CONTINE’ 

CALL GCMDV(20,7) 

WRITE (6,905) LIST 

CALL. GINKEY(I,MB) 

CALL GOLEAR(O,4) 

LIST = ' · 

CALL @СМОМ( 20,20) 

WRITE(5,805) LIST 

LIST = ‘THE DATA RECORDED FOR THIS EXPERIMENT EITHER:’ 
СА. асм0у(20,19) 

WRITE(6,805) LIST 

LIST = ‘ 

CALL @СМОМ(6,18) 

WRITE(5,805) LIST 

LIST = '1. REFLECTS THAT THE CHONGES IN CLERENT' 

CALL. GCHEOV (20, 17) 

ИНТТЕ (6,805) (ТЕТ 

LIST = 'MEASUREMENT DESCRIBED EARLIER DID NOT ОСОЛ: 
CALL ОСМОМ (20, 16) 

WRITE (6,805) LIST 

LIST = ‘AT A POINT(S), RELATIVE TO ECOOR, THAT WOLD’ 
СА ӘСМА (20, 15) 

WRITE(5,805) LIST 

LIST = ‘GENERATE ENOUGH DATA POINTS FOR THE USE OF THE’ 
CALL. GCHOV (20,14) 

WRITE(6,805) LIST 

LIST = ‘CENTRAL DIFFERENCE METHOD.’ 

CALL GCHDV (20,13) 

WRITE(6,805) LIST 
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LIST = 48: 
CAL ОСМОМ (59,11) 

WRITE(6,805) LIST 

LIST = '2. THE MODEL 351 DID NOT RECORD THE APPLIED’ 
CAL GN (20,9) 

WRITE(6,805) LIST 

LIST = ‘POTENTIAL PROPERLY. THIS CAN BE CORRECTED’ 
са. ОСМОМ (20,8) 

МВІТЕ (5,805) LIST 

LIST = ‘IF THE PROBLEM IS NOT CORRECTED THEN THE’ 
CALL ОСМОМ (20, 7) 

WRITE(6,805) LIST 

LIST = ‘CENTRAL DIFFERNCE METHOD WILL NOT WORK, AND THE’ 

CALL GCMOV(20,64) 

WRITE (6,805) LIST 

LIST = ‘CUBIC SPLINE METHOD MAY NOT WORK’ 

CAL GCOWXW(20,5) 

WRITE(6,005) LIST 

LIST = ‘PRESS ANY KEY TO CONTINE’ 

CAL ОСМОМ(20 5) 

WRITE(5,805) LIST 

CALL GINKEY(I,MB) 

CALL GOLEAR(O,4) 

LIST = ‘IF AN ATTEMPT IS MADE TO CORRECT THE’ 

CALL ОСМОМ (20, 17) 

WRITE(6,805) LIST 

LIST = ‘PROBLEM TE CORRECT METHOD IS DETERMINED’ 

CALL GOMOV(20, 16) 

WRITE (6,905) LIST 

LIST = ‘IF THE ATTEMPT IS UNSUCCESSFUL THEN THE’ 

CALL ОСМОМ (20,15) 

WRITE(5,805) LIST 

LIST = ‘CUBIC SPLINE METHOD SHOLLD EE USED. ' 

САД ОСМОМ (20, 14) 

WRITE(6,805) LIST 

LIST = ‘IF THE ATTEMPT IS SUCCESSFUL EITHER’ 

CAL ОСМОМ(20,15) 

WRITE(6,805) LIST 

LIST = ‘METHOD My BE USED.’ 

CALL GOMOV(20, 12) 

WRITE (6,805) LIST 

LIST = ‘IF THIS SET OF MESSAGES APPEPRS AGAIN’ 
CALL GOMDV(20,10) 

WRITE (4,805) LIST 

LIST = '@FTER ATTEMPTING TO CORRECT THE PROBLEM’ 
CAL QOXW(20,9) 

WRITE(6,805) LIST 

LIST = THEN TE CURRENT DENSITY VALLES ARE SUCH’ 
CALL ОСМОМ(20,8) 

WRITE (6,805) LIST 

LIST = ‘THAT NEITHER METHOD CAN BE USED’ 
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CALL ОСМОУ(20,7) 
WRITE(6,805) LIST 
LIST = 'PRESS ANY KEY TO CONTINE’ 
CAL ОСМОМ (20,5) 
WRITE(6,805) LIST 
CALL QGINKEY(I,MB) 
CALL GOLEAR(O,2) 
LIST = ‘1. DO NOT CORRECT THE PROHLEM. ' 
CALL GOMOV(25,13) 
WRITE(6,805) LIST 
LIST = ‘2. TRY AND CORRECT THE PROBLEM. ' 
СА ОСМОМ (25, 11) 
WRITE(6,805) LIST 
LIST = ЈЕМЕНА 1 СВА 2: 
CALL. GCMOV (32,9) 
WRITE(6,805) LIST 
CAL GCMOV(S9,8) 
READ(4,71°,ERR = 450, ТОЗТАТ = 34) М 
IF ((M .Gt. 2) „ОВ. (44 .NE. O)) THEN 
CALL OOCLEAR(O,4) 
LIST = · · 
CAL ОСМОМ (21,16) 
WRITE(6,805) LIST 
LIST = ‘YOU DID NOT ENTER AN INTEGER ОЕ’ 
CALL GOMOV(20, 15) 
МЕТТЕ( 6,905) LIST 
LIST = ‘VALLE 1 СК 2: 
CALL 9СМОМ( 57 , 14) 
WRITE(6,805) LIST 
LIST = ‘PRESS ANY KEY TO CUNTIN E: 
CALL ОСМОМ (25, 12) 
WRITE(6,805) LIST 
CALL ОСМОМ (59, 11) 
CALL QGINKEY(I,MB) 
GOTO 420 
ENDIF 
IF (M .ED. 1) THEN 
М = 1 
GOTO 1000 
ENDIF 


AN ATTEMPT IS BEING MADE TO CORRECT THE POTERTGIL VALUES SUCH 
THAT THEY ACCURATELY REFLECT THE EXPERIMENTAL POTENTIALS. 
BLOCKS SOO AND 600 CORRECT THE PROBLEM, BUT ONLY ONE (CF THE 
TWD ROCKS AE USED FOR ANY GIVEN FILE. 


IF (M .ED. 2) THEN 
DO SOO I = 1,N,2 
Y1 = MDATA(I,1) 
Y2 = MDATA(I+1,1) 
IF (Yi .ЕП. У2) THEN 


805 
810 


MDATA(I,1) = MDATA(I,1) – 0.000500 
№ =1 
ENDIF 
CONTINE 
PO ¿OO I = 2,N,2 
Yl = MDATA(I,1) 
у2 = мрата(1+1,1) 
IF (Yi „ЕО. Y2) THEN 
MDATA(I,1) = MDATA(I,1) — O.OO00SDO 
М = 1 
ENDIF 
CONT INUE 
IF (N9 .ED. 1) THEN 
CALL GOLEAR(O,2) 
LIST =" ° 
CALL GOMOV(30,14) 
WRITE (6,805) LIST 
LIST = ‘CORRECTION SUCCESSFUL’ 
CALL GOMOV(30,13) 
WRITE(5,805) LIST 
LIST = ‘USE EITHER METHOD’ 
CALL GoMov(20,11) 
WRITE (6,805) LIST 
LIST = ‘PRESS ANY KEY TO CONTINE’ 
CAL GOMOV(30,9) 
WRITE(6,805) LIST 
CALL. GINKEY(I,MB) 
GOTO 10 
ELSE 
CALL GOLEAR(O,4) 
LIST = ‘CORRECTION UNSUCCESSFUL ’ 
CALL. GCMOV (30, 13) 
WRITE(6,805) LIST 
LIST = ‘USE EITHER QUBIC SPLINE METHOD’ 
CALL GOMOV(30,11) 
WRITE(5,805) LIST 
LIST = ‘PRESS ANY KEY TO CONTINE’ 
CALL. GCMOV (30,9) 
АБІТЕ( 5,805) LIST 
CALL. QINKEY(I,MB) 
GOTO 1000 
ENDIF 


ENDIF 


FURMAT (ASO, \) 
FORMAT (I1) 
1000 END 
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XDAT : HE X VALLES OF THE FIRST GRAPH 

YDAT =: THE Y VALUES (F THE FIRST GRAPH 

XDAT1 =: THE X VALUES OG THE SECOND GRAPH 

YDATi1 : THE Y VALUES OF THE SECOND GRAPH 

TX > XDAT AS ORIGINAL PASSED AND SCALED. USED TO ALW 
THE USER TO RESCALE THE AXIS AND MAINTAIN ALL ПЕ 
VALLES ORIGINALLY PASSED 

TY : USED AS TX BUT FOR YDAT 

TAL > USED AS TX BUT FOR XDAT1 

ТҮ1 : ЕО А TX BUT FUR YDATI 


IN THIS SUBROUTINE ALL PROPIS PAE REPLIED TO BY THE LEE OF 
SFT KEYS. THIS MEANS THAT A CARRIAGE RETURN IS NOT USeD. 


SUBROUTINE GRAPH1 (CXDAT , YDAT ,N3,J „МАМЕ „ХОАТ? „УГАТ1 „МА, УМАМЕ „ППЕЗ, 


*51262. 01; ED SBETLES) 


REAL £4 XDAT (2000) , YDAT ( 2000) , XDAT1( 2000) , YDAT1 (2000) ,F ,F1,F2 
REAQLXA TX(2000) , TY (2000) ,TX1(2000) , TY1(2000) ,X2(1) , Y2(1) 
INTEBER32 K2,K3,K5,K5,L ,XL,XR, YL , YR, MB,01 ,62,C1,C2,M, IA(BO) 
INTEGER*2 1,J,N,K7,KB,K9, JCOL (10) ,JROW( 10) ,NL,NB,NS,NA ,N5,N6,N7 
INIEGER$2 L1,L2,L3,L4 

О-НЕҒҰЯСТЕН МЕМЕЖВ,ҮМЕМЕЯ5,ҒАСТА»,ХМОМЕЖ12,Т171Е%15,ТІТІЕ?2Ж52 
СНАРАСТЕВЖ2В, (АҒ5, ТІТІЕЗА5О 


N7 = 1 
M = O 

N = MAX(NS,N4) 
№ = NS — 3 
М = м – 5 


IF (N4 ЛЕ. O) N4 = 1 


THE ARRAYS ARE SCALED SUCH THAT THE Y AXIS VALLES ARE DISCERNIBLE 
FROM EACH OTHER 


F = MIN(ABS(YDAT(1)) ,ABSCYDAT(NS) ) ) 
Fl = MIN(ABS(YDAT1(1)) ,ABS(YDATI(N4) )) 
IF (N4 .GT. 2) THEN 
F2 = MIN(F1,F) 
El Se 
F2=F 
ENDIF 
IF (F2 .GE. 1.) THEN 
F2 = 1.0 
FACT = ‘х 1: 
ENDIF 
IF ((F2 .GE. .1) .ND. (F2 ЛТ L) ты! 
F2 — 1O 
РАСТ = “х 10” 
ENDIF 
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Н (ЕЛЕ БЕ Өг) AND. (F2 .LT. .L)) THEN 


F2 = 100. 
ҒЕСТ - “х 100 
ENDIF 


IF ((F2 „СЕ. .001) .AND. (F2 .LT. .01)) THEN 


F2 = 1000 

FACT = ‘х 1000’ 
ELSE 

IF (F2 .LE. 0.0010) THEN 

F2 = 1000 

РАСТ = "x 1000 

ENDIF 
ENDIF 
DO 6 I = 1,2000 

YDAT(1) = F24YDAT(I) 

ҮГРӨГ1(1) - Е2ЖҮРАГ1(1) 


CONTINUE 
= N3 
м = м 


THE PERMANENT ARRAYS ARE FILLED 


DO 10 I = 1,2000 


TX(1) = XDAT(1) 

ТҮ(1) - ҮГЕТ(І) 

ТХ1(1І) - XDAT1(1) 

ТҮ1(1) = YDATI(1) 
CONTINUE 


THE GRAPH AXIS AND WINDOW ARE DETERMINED 


XMAX = -1000 
XMIN = 1000 
YMAX = XMAX 
YMIN = XMIN 
ХМАХ1 = XMAX 
ХМІМІ - ХМІМ 
УМАХ1 = ХМАХ 
ҮМІМІ - ХМІМ 
хал - SO 
ЈСО 2 = 610 
JROW1 = 35 
JROW2 - ЗОО 


DO 200 I = 1,М5 

УМАХ = MAX(YDAT(I),YMAX) 
XMAX — MAX(XDAT (I) ,XMAX) 
XMIN — MIN(XDAT( I) ,XMIN) 
YMIN = MIN(YDAT(1),YMIN) 
CONTINE 

DO 210 I = 1,М4 

XMAX1 = MAX(XDAT1(I),XMAX1) 
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УМАХ1 = МӨХСҮПАТ16(1) ,ҮМӨХ1) 
XMIN1 = MIN(XDAT1(I),XMINL) 
YMIN1 = MIN(YDATL(1),YMINL) 
CONT INLE 

YFIN - (MAX(YMAX,YMAX1)) * .1 
YST (MIN(YMIN, YMIN1)) — .1 
XST (MIN(XMIN,XMIN1)) —. 1 
XFIN = (МАХ(ХМАХ,ХМАХ1)) + .1 
YMIN = YST 

YMAX = YFIN 

ХМАХ = XFIN 

XMIN - XST 

YORG - ҮМІМ 

XORG - XMIN 

XNAME = ‘LOG I A/cm2’ 

TITLE = *POTENTIODYNAMIC’ 
IDPT = O 

ҮСМЕНХ - 0.5 

ASPECT = 1.5 


|| 


l H H M 


|| 


THE GRAPH AXIS ARE DRAWN AND LABELED. TITLES ARE PLACED 


CALL СОМСЕЕ (16) 
CALL GPLOT(J00L1,JC0OL2, JROW1, JROW2, XMIN, XMAX , YMIN, YMAX, XORG, YORG. 


+ТОРТ , YOVERX , ASPECT) 


XMAJOR = (XFIN-XST)/10. 
XMINOR = O 

LABEL = -1 

МЕС = 2 

YMAJOR -(YFIN-YST)/10. 

YMINOR = O 

= 5 

= O 

IKOLOR = 1 

ІБҮМЕСІ --2 

KLRSYM-1 

QGCOLCOR(O,O) 

QSETUP (NDOTS, IKOLOR, ISYMBOL,KLRSYM) 
QXAXIS(XST, XFIN, XMAJCR , XMINCR , LABEL ,NDEC ) 
GYAXIS(YST , YF IN, YHAJCR , YMINCR, LABEL ,NDEC ) 
GXTICS(NKNT ,JCUL., JROMW) 

ОСТХТ (12, ХМАМЕ „2,3540,1,0) 
GYTICS(NKNT ,JCUL. , JROMW) 

ОБТХТ (5 „УММАМЕ „2,1,246,-1) 
QGTXT(6,FACT,2,1,162,—-1) 

QGTXT (15, TITLE, 2,254,335,0) 

QGTXT (8 NAME, 2,382,335,0) 

QSTXT (32, TITLE2, 2, 254,321,0) 

QGTXT (SO, TITLES, 2, 182, 307,0) 

= 2 

QSETUP (NDOTS, IKOLOR, ISYMBOL ,KLRSYM) 


$ 


 BREEREEERERRRE 
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CAL. GRTOI(XST,YST,L1,L2) 
CALL. GRTOI (XE IN, YF IN,L S, LA) 
CALL GLINE(L1,L2,L1,L4,2) 
CALL GLINE(L1,L4,L3,L4,2) 
CALL GLINE(L3,L4,L3,L2,2) 
CALL GLINE(L3,L2,L1,L2,2) 


THRE FIRST CURVE IS DRAWN. THE LINE COLOR IS CHANGED AND TIE 
SECOND CURVE IS DRAWN 


CALL GTABL(1,N5,XDAT,YDAT) 

КО ОК = 1 

CALL GSE TUP(NDOTS, IKOLOR, ISYMBOL,KLRSYM) 
CAI. GTOBL (1,N5, XDAT1, YDAT1) 

CALL (BEEP 

M=M+1 


IF TRE CURVES ARE OG DERIVATIVES AND A LINEAR REGION HAS ALREADY 
BEEN SELECTED, THE USER MAY DECIDE TO GENERATE ME TAFEL CONSTANTS 
BASED ON THIS PORTION OF THE OARE. 


ТЕ (СМ .GT. 1) .AND. (NA „БЕ.О)) THEN 
QUES = ‘GENERATE TAFEL CONSTANT’ 
CALL OGTXT(28,QLES,3,1,335,0) 


(RES = ' Y R N 72' 
CALL. GGTXT(28,0L ES, 5,1,.520,0) 
CAL GBEEP 


CALL. GINKEY(I,MB) 

IF ((MB .EQ. 89) „ОК. (МВ „ЕО. 121)) THEN 
Бїс эн - 

CALL OGTXT(25,QLES,3,1,335,0) 

CALL GGTXT(28,0 E5,5,1,320,0) 

QES = ‘IS THIS TE ANDIC CR' 

CALL GGTXT(28,10 E5,3,1,335,0) 

QUES = ‘CATHODIC BRANCH’ 

CALL GGTXT(28,0 £5,3,1,320,0) 

аЕ5 -“ АСС 2: 

CALL GGTXT(28,QLES,3,1,305,0) 

CALL. GEEEP 

CAL QINKEY(I,MB) 

IF ((MB .EQ. 65) .CR. (MB .EQ.97)) THEN 
Аі = кв 


IF ((MB .EG. 67) .ОН. (MB .ЕО. 99)) TEN 
Ci = KG 
C2 = K7 
GOTO 390 


я 


GOTO 381 
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390 


410 


420 


ENDIF 

ENDIF 

ІР ССМ .ЕП. 73) Ок. (МВ -ЕС ГЭ 002 85 
СОТО 590 

965 = 
GOTO 580 

ENDIF 

ENDIF 


BY THE USE OF THE CURSOR KEYS ON THE NUPERIC KEYPAD THE USER 
SELECTS A NEW SET OF AXIS TO BE PLOTTED 


ЕТ 
CAL ОБТХТ(28,2Е5,5,1,355,0) 
CAL GGTXT(2B,Q ES,3,1,320,0) 
CALL GGTXT(2B,QL ES, 3,1 ,:305,0) 
QUES = ‘NEW AXIS Y RN? 
CALL OGTXT(23,QLES,3,1,335,0) 
CALL СЕЕЕР 
CALI. GINKEY ( I,MB) 
IF ((MB .EG. 89) .CR. (МВ .ЕО. 121)) THEN 
GOTO 400 
ELSE 
IF ((MB .ЕП. 78) .CR. (MB .EB. 110)) THEN 
GOTO 800 
ELSE 
GOTO 390 
ENDIF 
ENDIF 
QUES = ° 
CALL OGTXT(28,Q@ES,1,1,335,0) 
CALL QGTXT(28,QLES,1,1,320,0) 
QES = ‘NEW,LAST,OR ORIGINAL (N,L,0)’ 
CALL ОБТХТ (28,01Е5,5,1,555,0) 
CAL СВЕЕР 
CALL. GINKEY(I,MB8) 
ЇЕ ((MB .EG. 79) .CR. (МВ .ЕО. 111)) THEN 
DO 410 І - 1,N3 


XDAT(I) — TX(I) 
УрАТ(1) = ТУ(1) 
CONTINE 


DO 420 I = 1,N4 
храта (1) = TX1(1) 
YDAT1(1) = TY1(I) 
CONTINUE 
M = O 
GOTO 5 
ENDIF 
IF ((MB .EG. 76) .CR. (MB .EG. 108)) THEN 
XST = XST1 
XFIN = XFIN1 
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YST = YST1 
YF IN — YFINI 
NS = MS 
№ = М6 
GOTO 530 
ENDIF 
IF (TE TEOT) -UR 
XST1 = XST 
YST1 = YST 
XFIN1 - XFIN 
YFIN1 = YFIN 
№ = № 
мМ = № 
(AES = 
CALL 


(MB 


„Ей. 


QGTXT(28,Q01E£5,.5,1,. 535,0) 


110)) THEN 


(AES = °TO THE LOWER LEFT HAND’ 


Y2(1) 


GGTXT(ZB,QLES, 3,1 ,:320,0) 


“CORNER, HOME TO EXIT’ 


„ЕО. 72) THEN 
GEMOVR (0,5) 
900 


ЕО. 75) TÆN 
Н МОУНС-5,0) 
СОО 


-EQ. 77) THEN 
СА МОМЯ(5,0) 
900 


-EQ. 80) THEN 
GLMDVR(0,-5) 
900 


-EQ. 71) THEN 
G_LREAD(XL, YL) 


GGTXT (28,GLES,3,1,305,0) 
ОВЕЕР 

СА 5ТАТ(1) 

G TYPE(3,3) 

G MV (320,175) 
GINKEY(I,L) 


GL MOVR(S20-XL ,175-YL) 
QITOR(XL,YL,XST,YST) 
X2(1) = XST 


= YST 


CALL GSETUP(5,3,43,3) 
CALL ОТАН (0,М7,Х2,Ү2) 
СОТО 510 


ELSE 


GOTO 430 
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ENDIF 
510 QUES = ‘USE THE (GRSUR KEY ТО МХБ” 

CAL GGTXT(28,Q ES,3,1,535,0) 

QUES = ‘TO THE UPPER RIGHT HAND’ 

CALL QGTXT(23,QES,3,1,320,0) 

GLUES = ‘CORNER, HOE TO EXIT’ 

CALL GGTXT(28,QUES,3,1,305,0) 

CALL GEEEP 
520 CALL. QINKEY(I,L) 

IF (L .EQ. 72) THEN 

CAL G MNR(O,5) 


БОТО 520 

ENDIF 

IF (L .EQ. 75) THEN 
САЛТ С МОМЯ(-5,0) 
GOTO 520 
ENDIF 

IF (L .EU. 77) THEN 
CALL СА МОМҒ(5,0) 
GOTO 520 
ENDIF 

IF (L .ED. 80) TFEN 
CALL GLOVR(CO, 75) 
GOTO 320 

ENDIF 

IF (L .EQ. 71) THEN 
CALL GLREAD(AR, YR) 


CALL. GL MOVRCSZO-XR, 175- YR) 
CALL. OITOR(XR, YR, XF IN, YF IN) 
X2(1) = XFIN 
Y2(1) = YFIN 
CALL атан (0,М7,Х2,Ү2) 
CAL G STAT(O) 
GOTO 530 
ELSE 
GOTO 510 
ENDIF 
ENDIF 
530 IF (YNAME .ЕО. ‘VOLTS’) THEN 
GOTO 780 
ENDIF 
Kb 
K7 
DO 600 I 
X = TX(I) 
Y-TY(I) 
IF (((X „СТ. ХТ) „ИМО. (X .LT. XFIN)) .eND. ((Y .GT. YST) 
+ .AND. (Y .LT. YFIN))) THEN 
Ко I 
GOTO 505 
ENDIF 


1 
NS 
= 1,N3 


6005 


670 


700 
715 


770 


CONTINE 
DO 610 I = K6, N3 
X = TX(I) 
Y — TY(I) 
IF (((X .LT. XST) .OR. (X .GT. XFIN)) .OR. ((Y .LT. YST) 
.CR. (Y .GT. YFIN))) THEN 
K7 = I 
GOTO &30 
ENDIF 
CONTINUE 
№ = K7-K& 
00 670 I = K&,K7 
YDAT(I-Ké*1) = TY(1) 
XDAT (I-KS+1) = TX(I) 


CONTINE 


|| 
P 


80036 
ШІ 


IF (((X „СТ. XST) -AND. (X .LT. XFIN)) .AND. ((Y .GT. YST) 
„РМО. (Y .LT. YFIN))) THEN 
K8 = I 
БОТО 715 
ENDIF 
CONTINE 
DO 720 I = K8,N4 
X = ТХ1(1) 
Y = ТҮ1(1) 
IF (((X .LT. XST) .CR. (X .GT. XFIN)) .CR. ((Y .LT. YST) 
.CR. (Y .GT. YFIN))) THEN 


K? = I 
БОТО 730 
ENDIF 
CONTINE 


DO 770 I = KB,K9 
YDAT1(I—K8+1) = TY1(1) 
XDATi(I-KB*-1) = TX1(1) 

CONTINUE 

Nb = K? — KB 

GOTO 225 

N5-0 

N5-0 

DO 785 I = 1,N3 

Х = ТХ(1) 
У= ТУ (1) 
IF (((X .GT. XST) .6ND. (X .LT. XFIN)) .6ND. ((Y „ВТ. YST) 
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* -AND. (Y ЦТ. YFIN))) THEN 
№ = № +1 
XDAT(N5) = TX(T) 
YDAT(NS) = ГИГ) 


ENDIF 
735 СОВА МЕ. 
DO 790 I = 1,М4 
X = TX1(I1) 
Y-TY1(I) 
IF (((X .GT. XST) .^ND. (X .LT. XFIN)) .&ND. ((Y .GT. YST) 
* .CND. (Y .LT. YFIN))) THEN 
№ = № +1 
XDATI(NS) = ТХЇ(1) 
ҮРАТ1 (Мб) - ТУ1(1) 
ENDIF 
790 CONTINUE 
M = O 
GOTO 225 


800 GLUES = ‘DO YOU WANT TO PLOT’ 
CALL GGTXT(23,QLES,3,1,335,0) 
QUES = ‘THIS GRAPH Y CR N ?' 
CAL OGTXT(28,QES,3,1,320,0) 
CALL СВЕЕР 
CAL OINKEY(I,MB) 
IF ((MB .EB. 89) .OR. (MB .EQ. 109)) THEN 
GOTO 830 
ELSE 
IF ((MB .EG. 78) „ОВ. (МВ -ЕО. 110)) THEN 
вото 850 
El SE 
GOTO 800 
ENDIF 
ENDIF 
830 CALL PLOT(XDAT,YDAT,NS, XDAT1L, YDATL,NS, XNOME, YNEME, TITLE, 
+XMAX, XMIN, YMAX, YMIN, XST, XFIN, YST, YFIN, XMAJOR, YMAJOR, TITLE2, 
FACT, NAE) 
IF (YNAME .ЕП. 'MLTS') THEN 
COTO 390 
ENDIF 
ЈЕ (УМЕ „.ЕО.  'Dv/Di') THEN 
СОТО 380 
ENDIF 
850 DO 860 I = 1,2000 
ТХ(1)-0 
ТҮ(1)-0 
TX1(I) = O 
ТҮ1(1) = О 
XDAT(I) = O 
XDAT1(I) = О 
YDAT(I) = 0 
Үрат1(1) = о 
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CONTINUE 

CALI. GOIODE ( 16) 
САЦ. СЕМООЕ( 3) 
FORMAT CÓ) 

END 
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THIS SUBROUTINE PLOTS THE GRAPHS AS SHOWN ON TRE SCREEN DURING 
TRE CALLING SUBROUTINE GRAFH. 

THE MAJORITY OF THE PROGRAM IS: INVOLVED IN SETTING UP ARGLEMENTS 
ТО HE PASSED TO SUBROUTINE CONTAINED WITHIN ТРЕ PLOTMATICS 
LIBRARY. THIS MANUAL SHOULD BE CONSULTED IN ORDER TO DETEMINE 
THE USE OF THE ARGLEMENTS OTHER THAN THOSE. DEFINED BELOW 


ADOT : THE ARRAY CONTAINING THE X VALLES OF THE FIRST CURVE 
TO BE PLOTTED. 


YDAT = THE FIRST CURVE’S Y VALES 

XDAT1 : THE ARRAY CONATINING THE Y VALUES OF THE SECOND ORVE 
ТО ВЕ PLOTTED. 

YDAT1 : THE SECOND CURVE’'S Y VALES 

TIME = TFE TITLE Te Peer 

7" ПЕ2: AN DISCRIBTICN OF THE GRAPH 

ANC" : THE X AXIS TITLE 

SAPPLE: THE Y AXIS TITLE PROPERLY POSITION AND CATENTATED 


SUBROUTINE PLOT(XDAT, YDAT,NPT1, XDAT1, YDAT1,NPT2, XNAME , YNAME 
+, TITLE, XMAX, XMIN, YMAX, YMIN, XST, XFIN, YST, YFIN, XMAJOR, YMAJOR 
+, TITLE2,FACT, NAME) 

REALX4 XDAT (2000) , YDAT ( 2000) , XDAT1(2000) , YDAT1 (2000) 

INTEGERX2 JCOL(10),JROW(10) ,NPT1,NPT2 

INTEGER 15,14,15,16,17,18,111,112,M8 

CHARACTERXB МАМЕ , YNCIE XS, XNCMEX 11, TTTLEX15, 0E E12, TT TLE2X 2 
+, FACT#6,04%1, GUES*29 


да = 
SAMPLE = YNAME//A4//FACT 
ХОКб = хот 
YORG = YST 


DTX = ABS(XMAX-XMIN) 

DX = DTX /MAX(NPT1,NPT2) 
DTY = ABS(YMAX-YMIN) 
ТОРТ = © 

NKNT — 5 

NDOTS — O 


ISIZE - 
ІСТЕР - 
ТРОКТ = 
CALL ZINITCISIZE, ISTEP, IPORT, BATCH) 
JCL = 20 

JCR = 620 

ICL 10 

ICR 320 

15 =. 12 

I4 - 150 


| 
= N = 
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ГЭ 
16 
I7 
I8 


Hg 


111 - 


222 


I12 = 182 


CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
ХМІМІ 
ХМАХ 1 
УМАХ1 
ҮМІМІ 


PERRA 


: 


23431131433331314142134343343: 


сами (аа „та. „16) 

ZCONV (JOR, [OR, 14) 

700МУ (13, 14,16) 

ZOONV( 15, 16,14) 

ТООМУ( 17, 18, 16) 

ТСОММ(111, 112,16) 

-XMIN —.2XDTX 

= XMAX+ .1*(DTX) 

YMAX + .1XDTY 

YMIN — .1XDTY 

= 2 

ZLOT(JCL,JCR, ID I ICR, XMINI ,XMAX1 , YMINI , YMAX1 , XORG , YORG) 
ZRTOI (XST , YST ,L.1,L 2) 

ZRTOI ( XFIN, YF IN,L 5,L 4) 

15ЕПР(-1,1,0,1) 

= 1 

ZXAXIS(XST, XFIN, XMAJOR ,MINOR,LABEL ,NDEC) 

ZYAXIS(YST , YF IN, YMAJOR , MINOR, LABEL „МОЕС) 
ZXTICS(10,JCO_,JROW) 
ZYTICS(10,JCOL,JROW) 
ZLINE(L1,L2,L1,L4,1) 
ZLINE(L1,L4,L3,L 4,1) 
ZLINE(LS,LA4,L3,L 2,1) 
ZLINE(L3,L2,L1,L 2,1) 
ZPU 
QGTXT(28,QLES,3,1,335,0) 
ӨЗТХТ(28,(ПЕЗб,3,1,3520,0) 
ӨЗТХТ(28,ЕЗ,3,1,305,0) 

= ‘IF DESIRED INSERT A NEW PEN’ 
QGTXT (28, QLES,3,1,335,0) 

= ‘THEN PRESS ANY KEY’ 

QGTXT (28,GLES,3,1,320,0) 
ОВЕЕР 

QINKEY(I,MB) 

ZTABL (1,NPT1 , XDAT , YDAT) 

ZPU 

= ‘IF DESIRED INSERT ñ NEW PEN 
OGTXT (28,QG.ES,3,1,335,0) 

= ‘THEN PRESS ANY KEY’ 

QGTXT (28,QES,3,1,320,0) 
GEEEP 

GINKEY (I, B) 

ZTAOEL. (1,NPT2, XDAT1, YDAT1) 
ZPU 
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АТ 


р 
Т 


4343431341333: 


GINKEY(I,MB) 
= ‘iF DESIRED INSERT A NEW РЕМ' 
ОБТХТ (28,01Е5,3,1,335,0) 
= ‘THEN PRESS ANY KEY’ 
QGTXT (28,GLES,3,1,320,0) 
СВЕЕР 

гетхта(11,хчаме ‚ 1) 
ZSETUP(2,1,0,1) 
Z1W(1,1,1750,1400,0) 
ZPTXT(15,TITLE, 1, I5, I6) 
ZPTXT (8, NOME, 1,111,165) 
ZPTXT (32, TITLE2, 1,15, 18) 
ТЕХТХТ(О,0,0,0,0,1) 

ZPU 

гра(15,14) 

ZLABEL (12,SAMPLE) 

ZFINIS 
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THIS IS THE FIRST SUBROUTINE CALLED WHEN USING IHE CLBIC 
SPLINE METHOD. 


THE ARRAY CONTAINING POTENTIAL AND LOG CURRENT DENSITY 

THE ARRAY CONTAINING THE LOG CURRENT VALLES FOR THE 

CATHODIC BRANCH 

YDAT =: THE ARRAY CONTAINING THE DERIVATIVES FOR THE CATHODIC 
BRANCH 

XDAT1 : SAME A XDAT BUT FOR THE ANODIC BRANCH 

YDAT1 =: SAME AS YDAT BUT FOR THE ANODIC BRANCH 

DIV : ПЕ ARRAY CONTAINING THE DERIVATIVES WHICH IS BEING 

RETURN FROM SUBSEGLENT CALLED SUBROUTINES 


MDATA 
XDAT 


THE INTEGER VARIABLES ACTUALLY USED IN THIS SUBROUTINE CHE 
USED TO DEFINE THE REGIONS WITHIN BOTH BRANCFES AVAILABLE 
FOR USE 


SUBROUTINE G_OPE1 (MDATA,N,J,NAME,A1,A2,C1,C2) 

REAL X4 XDAT(2000) , YDAT ( 2000) , XDAT1(2000) , YDAT1(2000) 
REAL XB MDATA(2000,2) ,DIV(1000) 

INTEGER32 I,J,N,K,MI ,NS,NA,N1,K1,01,02,C1,C2 
CHORACTER NOEXB,YNOPE$5,TITLE2X32, TITLE SK5O 

Ni = № 

YNAME = “Ім/Пі” 

TITLE? = 'TCFEL SLOPES: CUBIC SPLINE’ 

TITLE3=" 


THE CATHODIC BRANCH IS CHECKED FOR AN INCREASE IN LOG 
CURRENT DENSITY WHICH CALD САБЕ THE CAOLATION CF ñ 
DERIVATIVE OF INFINTE SLOPE 


DO 2O I = J-1,2,-1 
IF (MDATA(I-1,2) -LE. FDATA(1,2)) GOTO 25 

CONTINUE 

GOTO 26 

K = I 

MI = J-9 

J1 = K 

МЗ ABS (K-MI ) 

М = 7 

CALL CSPLIN(MDATA,N,M,DIV,K,NS) 

DO 100 I = 1,N3-3 

XDAT (I) = моата(1-4-1,2) 

УПАТ( 1) = DIV(N3-I+1) 

CONTINE 
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126 


TE ANDIC BRANCH 15 CHECKED FOR LOG CURRENT VALUES WHICH 
COULD CAUSE THE CALOCULATION OF A CERIVATIVE GF INFINITE SLOPE 


DO 120 I = J+1,Ni 
IF (MDATA(I,2) -LE. MDATA(I-1,2)) GOTO 125 


CONTINUE 

K = MIN((J*290) , (N1)) 
GOTO 126 

K = I 

J1 = J 

MI = K — 9 

№ = PPS(MI—J) 


CALL CSPLIN(MDATA,N,M,DIV,J,N4) 
DO 200 I = 1,N4 


XDAT1L(1) = MDATA(J1+I+4,2) 

YDATL(I) = ОКТ) 

CONTINE 

CALL. GROPH1 (XDAT , YDAT ,N5, J NOME , XDAT1 , YDAT1 ,NÀ , YNGYE , TTTLE2, 
+01 ,A2,€1,C2, TITLES) 

END 
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THIS SUBROUTINE SETS UP THE MATRIX CF CCEFFICIENTS ТО 
ВЕ SOLVED IN ORDER ТО DETERMINE THE DERIVATIVE USING THE 
CUBIC SPLINE METHOD. 


MDATA : THE ARRAY CONTAINING POTENTIAL AND LOG CURRENT DENSITY 
CMAT : THE 7x7 MATRIX CONTAINING THE COEFFICIENTS 
YMAT =: THE 7x1 VECTOR CONTAINING THE RIGHT HAND SIDE. OF 


THE SYSTEM GF EQUATIONS. WHEN RETURNED FROM THE 
SOLVING SUBROUTINE IT CONTAINS THE SOLUTION. 


GUBROUTINE CSPLINE (MDATA,J,M,DIV,K,MI ) 
REAL’8 XDAT(1000), YDAT( 1000) ,MDATA(2000,2) , YMAT(7) 
REAL*B CMAT(7,7),AI,BI,CI 
REZL $8 DIV(1000) 
ІМТЕСЕКЖ2 М,МІ,К,Ј 
ГО 101 = 1,7 

00 5 Ј = 1,7 

OMAT(I,J) = O.DO 


CONT TINUE 

CONT IN E. 

М- 7 

K = K-1 

DO 200 М = 1,MI 

к-к+1 

CMAT(1,1) = 2.D0O8(MDATA(K+2,2)-MDATA(K,2) ) 

CMAT (1,2) = MDATA(K+2,2) — MDATA(K+1,2) 

CMAT(7,7) = 2. D0 (MDATA(K*8,2) -MDATA(K*6,2) ) 

Смат (7,6) = MDATA(K+8,2)—-MDATA(K+7,2) 

YMAT(1)- { (МОАТА(К+2,1)-МОАТА(К+1,1))/(МОАТА(К+2,2) -МОАТА 
+(К+1,2))) — ((МОАТА(К+1,1)-МОАТА(К,1))/ (МОАТА(К+1,2)-МОАТА 
+(К,2))) 

YMAT(1) = YMAT(1)*6. 

ҮМАТ (7) = ( (MDATA(K+8, 1 )—-MDATA(K+7,1) )/ (MDATA(K+8, 2)-MDATA 


+(K+7,2))) — ((MDATACK+7,1)—-MDATA(K+6, 1) )/ (MDATA(K+7,2)—-MDATA 
*(K*6,2))) 
ҰМАТ(7) - ҮМАТ(7)ж6. 
DO 100 I - 2,6 
K = K +1 
CMAT(I,I) — 2.DOX( (MDATA(K*2,2) HMDATA (K,2) )) 
СМАТ(Т, 1+1) = МОАТА(К+2,2) МОАТА(К+1,2) 
CMAT(I,1-1) = MDATA(K+1,2) — MDATA(K,2) 
УМАТ(Т) = ((МОАТА(К+2,1)-МОАТА(К+1,1) ) /СМАТ(Т,1+1)) 
+— ( (МОАТА(К+1,1)МОАТА(К,1) ) /СМАТ(Т,1-1)) 
УМАТ(Т) = УМАТ( Г) же. 
CONTINE 
к-к-5 
CALL. LINSY1(CMAT , YMAT ,M) 
АІ - ҮМАТ(5) - ҮМАТ(4) 
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AI/((6.DO) *( (VDATA(K*4,2) H7DATA(K 35,2) ))) 
YMAT (4) /2.DO 
СІ (МОАТА(К+4 ,1) МОАТА(К+3,1))/ (МОАТА(К+4 ,2) МОАТА(К+$,2)}) 
CI CI-(1/6) *( ( (23YMAT (4) )+YMAT (5) ) 8 (MDATA(K+4 ,2)-MDATA(K+ 
+5,2))) 
DIV(N) = CI 
200 CONTIN E 
END 


AL 
BI 
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THIS SUBROUTINE SOLVES THE 7x/ SYSTEM OF ECAATIONS GENERATED 

BY CSPLIN.FOR. THE PROGRAM WAS ORIGINALLY FROM THE TEXTBOOK 
“FORTRAN 77 for Engineers and Scientists” BY NYHOFF AND 

LEFSTMA, PAGES 264-265. THE MAIN MODIFICATION IS THE ELIMINATION 
OF THE VARIABLE DIMENSIONING OF ARRAYS PRESENT IN THE ORIGINGL 


PROGRAM, 


SUBROUTINE LINSY1(COMAT, YMAT,N) 
INTEGER N,PIVOT 
REAL*B AUG(7,9), COMAT(7,7), YMAT(7),X(7) „МА Т,ТЕМР 
D0 20 I=1,N 
DO 10 J = 1,N 
А5(1,Ј) = 
CONTINUE 
CONTIM E 
DO 3O I = 1,N 
AUG(I,N+L) = YMAT(1) 
CONTINUE 
DO 70 I = 1,N-1 
IF (AUG(I,I) .EQ. O.) TEN 
PIVOT = O 
J = 1+1 
IF ((PIVOT .EQ. О) AND. (J .LE. N)) THEN 
IF (AUG(J,1) .МЕ. O.) PIVOT = J 
J = Ј+1 
GOTO 35 
ENDIF 
IF (PIVOT .EQ. O) THEN 
PRINT*, ‘SINGULAR’ 
STOP 
ELSE 
DO 40 J = 1,N+1 
ТЕМР = AUG(I,J) 
АШБ(1,Ј) = AUG(PIVOT,J) 
AUG(PIVOT,J) = TEMP 
CONTINUE 
ENDIF 
ENDIF 
DO 60 J = I+1,N 
MLT = +¥%G(J,1)/AUG(I,1) 
DD SO K = I,N+1 
АВ 9 ,„К) < АБК) + МАТЗАЮ(Т,К) 
CUNTIN E 
CONTINUE 
CONTINUE 
X(N) = AUG(N,N+1)/AUG(N,N) 
DO 90 J = N-1,1,-1 
X(J) = AUG(J,N+1) 


COMAT(1,J) 
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D 8O K - J*1,N 
X(J) — X(J)-(AJ)GC  J,K) &X(K)) 


X(J) 2 X(J)/AJG(J,J) 


CONTINUE 

DO 100 I = 1,N 
YMAT(I) = X(T) 
CONTINE 

END 
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THIS SUBROUTINE СОМЯНТ5 А МҒЕНІС УҮКІНЕ (Х6) CF FORMAT 

+/— XX.XXXX TO A CHARACTER VARIBLE WHOSE WIDTH IS 8 CHoROCTEPS. 
THIS IS NECESSARY AS THE PLOTMATICS LIBRARY IS NOT EQUIPPED WITH 
A SUBROUTINE FOR PLOTTING NVPERIC VARIBLES. 


SUBROUTINE CONV(A7,X&) 
INTEGER 01,02,03,04,05,06 

CHARACTER*1 (0:9), С1,С2,С5,С4,С5,С5,С7,СВ,А7%8 
66(0) = · 
ASé(1) 
66(2) 
^63) 
66(4) 
66(5) 
66(6) 


е 


ин 
чичсивзямно 


$ 
o 
` H H H 


(X6 „ЕТ. 0.0) С1 = *— 
ABS( X6) 

X6/10. 

IFIX(X1) 

- X6-(10.X01) 

IFIX(X2) 

X6 — (10.301) - 02 
X3/.1 

- IFIX(X3) 

= (Х3-05)/.1 + .001 
- ЇЕЇХ(Х4) 

- (Х4-04)/.1 

- ІБІХ(Х5) 

- (X5-05)/.1 

- ІБІХ(Х7) 

10 I = 0,9 

ТЕ (01 .ЕП. 1) C2 
IF (O02 .ED. I) СЗ 
IF (O3 .ED. I) CA 
ТЕ (04 „ЕД. 1) CS 
IF (05 .ЕП. 1) Сё 


>x О > >< 
N P = с T 
| HO H H 


SFIGASSAAGA 
И 


AS(T) 
AS(T) 
AS(T) 
AS(T) 
AGI) 
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IF (06 EQ. I) CZ = 66(1) 
CONTINE 
св 
А/ - С1//С2//С5//СВ//С4//С5//С6//С7 
RETURN 
END 
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THIS SUBROUTINE SETS SOME (F THE CHORACTER. STRINGS USED 
WHEN USING THE CENTRAL DIFFERENCE METHOD 


SUPRDOUTINE G OPE (MDATA,N, J NAME ,A1,AZ,C1,C2,N1,N2,N3,N4) 
FEAL *4 XDAT (2000) , YDAT (2000) , XDAT1 (2000) , YDAT1 (2000) 
КЕЯ 8 МОАТА( 2000,2) 
INIEGERX2 I,J,N,01,02,C1,C2 
СНАРАСТЕК NOG"EX8,YNOI'E$5,TITLE2*32, LIST€&SO,TITLE.SXSO 
YNAME = °Dv/Di’ 
TITLE? = TAFEL SLOPES: CENTRAL DIFFERENCE’ 
TITLE=’ 
CALL DATADEL (MDATA,J.YDAT,YDAT1, XDAT , XDAT1,N2,N1) 
CALL GRAPH (XDAT , YDAT ,NS,J, NAME, XDAT1, YDAT1,N4, YNAVE, TITLE?, 
+A1 ,A2,C1,C2, TITLES) 
END 
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x THIS GUBROUTINE CONDUCTS NUMERICAL DIFFERENTIATION USING 
THE 4 POINT CENTRAL DIFFERNCE ГЕТН. 


ж 


4 MDATA = ARRAY CONTAINING POTENTIAL AND LOG I CENSITY 
* я ОРЕ = ARRAY CONTAINING DE/D(LOGI) AND LOG I 
x ЕЯ_ОРЕ = CENTRAL DIFFERNCE DIFFERENTIATION 
3 ECORR1 = UPPERLIMIT CF ANODIC/CATHODIC DISTRIBUTION 
Ж ECURR2 = LOWERLIMIT OF ANDDIC/CATHODIC DISTRIBUTION 
x N = NUMBER OF DATA POINTS 
3 3 = LOCATION WITHIN MDATA OF ECORR 
4 I = CATHODIC BRANCH LOOP COUNTER 
x K = ANODIC BRANCH LOOP COUNTER 
Ж H = DIFFERENTIAL VALUE IN POTENTIAL 
4 11 = VALLE OF LOG I AT E—2H 
k I2 = VALLE OF LOG I AT E-H 
x 15 = VALLE OF LOG I ATE 
x 14 = VALLE OF LOG I AT E+H 
x I5 = VALLE OF LOG I AT Е+2н 
* DEDI = SLOPE VALE O 13 
SUBROUTINE DATADEL (MDATA,J,YDAT,YDAT1, XDAT,XDATL,N2,N1) 
REAL*8 MDATA( 2000, 2) ,X1,X2, X3, X4, X5,H,FSLOPE, DEDI 
REAL #4, XDAT1(1000), YDAT1(1000), YDAT(1000) , XDAT( 1000) 
INTEGER*2 J,K,I 
FSLOPE (Н,Х1,Х2,Х4,Х5)-1Х1-(8.008Х4)-1(8.008Х2)-ХЭ)/(12.ПО8Н) 
C PERFORMING CALCULCATIONS ON CATHODIC BRANCH 
K =O 
DO 100 I = J—2,N1+2,-1 
K = K+1 


H = MDATA(I,1) — MDATA(I-1,1) 
X мрата(1–2,2) 
MDATA(I-1,2) 
X3 = MDATA(I,2) 
MDATA( 1+1,2) 
MDATA(I*2,2) 
DEDI = 1.DO/FSLOPE(H, X1, X2,X4,X5) 
YDAT(K) = DEDI 
XDAT(K) = MDATA(I,2) 
100 CONTINE 
1-0 


> 
N = 
| | 


x 
Ол 
ІІ 


Е PERFORMING CALCULATIONS ON TMTE ANDDIC HRANCH 


200 DO 300 K = J+2, N2-2 
I 12-21 
H = MDATA(K+1,1) — MDATA(K,1) 


= 
— 


1.56 


X1 — MDATA(K-2.2) 
XZ = MDATA(K-1,2) 
х5 = мрата(к,2) 

ХА = MDATA(K+1,2) 
XS = MDATA(K+2,2) 


DEDI = 1.DO0O/FSLOPE(H,X1,X2,X4,X5) 
YDAT1(I) = DEDI 
XDAT1(I) —- MDATA(K,2) 
хо CONTIN E 
400 RE TURN 
END 
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